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SUMMARY 


Analytic methods are developed for calculating blade loads and shaft- 
transmitted vibratory forces in stiff bladed hingeless rotors operating at 
advance ratios from |*i = . 3 to |i = 2.0. 

Calculated shaft harmonic moments compared well with experimental 
values when the blade first flap frequency was in the region of two- 
per- revolution harmonic excitation. Calculated blade bending moment azi- 
muthal distributions due to changes in cyclic pitch agreed well with ex- 
periment at radial stations near the blade root at values of the ratio of 
first flap frequency to rotor rotation rate from 1.5 to 5 - 0 . At stations 
near the blade tip good agreement was only obtained at the higher values 
of frequency ratio. 

A compendium of experimental shaft transmitted force and blade loads 
data for two different rotor systems is included in Volume II. 
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INTRODUCTION 


Within the spectrum of VTOL aircraft types there is a class that must 
have a low downwash velocity for efficient hovering and to avoid excessive 
disturbance to the surface and persons standing below. It also must possess 
high speed, long range, and flight characteristics that are not fatiguing 
to the pilot or passengers even in gusty conditions. It must also possess 
a low initial, cost and inexpensive maintenance to be commercially or mili- 
tarily viable. 

The low downwash velocity in hovering flight demands a low disk loading 
or large disk area per unit weight. At present this requirement can be met 
only by large diameter rotors. The need for rotors designed for high for- 
ward speed, low drag, rapid highly damped response to control motions, mild 
response to gusts, and low vibration level has stimulated a search for a 
simpler, more inherently stable rotor system than found on most helicopters 
in operation today. 

One approach to meeting these requirements has been the development of 
the "hingeless" rotor. The hingeless rotor differs from those found in most 
currently flying helicopters in that its blades possess only bearings for 
feathering motion, which is controlled rather than free, and no flapping or 
lead-lag hinges. The blades of most conventional helicopters are essentially 
gimballed near the mast or they may be joined rigidly together and the whole 
rotor gimballed to the top of the mast (teetering rotors). 

Two main advantages of the hingeless or c ant ilever-b laded rotor are its 
ability to apply hub moments through the shaft to the body rapidly and of 
generally greater magnitude than available from articulated rotors, and its 
making available larger and more appropriately aligned damping forces due to 
body pitch and roll rotation rates. These potential advantages, of course, 
have long been appreciated in the helicopter world for the cantilever-bladed 
rotor. They could not be taken advantage of in the past, however, due to 
the extreme gust sensitivity of such a system and the extreme and variable 
cross-coupling of the cyclic and lift producing controls. 



With the development, over the past fifteen years, of the directly 
geared gyroscope-stabilized feathering system, gust sensitivity and controls 
cross-coupling have been greatly reduced and rapidly responding hingeless 
rotor helicopters with highly stable body modes have been developed and demon- 
trated. 

The development of the hingeless rotor gyroscope-stabilized system, how- 
ever, has not been without its problems, most of which stemmed from the lack of 
significant damping in the blade in-plane modes. As opposed to the flapping 
modes which are highly damped by aerodynamic forces, the in-plane aerodynamic 
damping is small. Under near 2P resonant conditions, the in-plane forces and 
motions can combine with the flapping deflections and forces to nullify the 
hub moment feedback path to the gyroscope and produce instability. This 
particular problem has been solved in two ways. First, the in-plane reson- 
ance has been kept to a low rotor rpm and tight control has been maintained 
on the rpm to avoid the dangerous region. The second solution has been to 
substitute direct hub moment feedback for the hub moment feedback via feath- 
ering moments, thus eliminating the in-plane motion sensitive part of the 
feedback path. 

Another in-plane resonance problem occurs in rotors with more than three 
blades. Pilot-induced oscillations can lead to high stresses in the reaction- 
less modes. These modes are particularly lightly damped since there is no 
shaft motion present in the mode and therefore no damping from the body. 

Again, strictly keeping away from the resonance has been the solution. 

One more possible consequence of light in-plane damping, which can be 
avoided by careful distribution of mass and stiffness on either side of the 
feathering hinge, is flap-lag instability. This is caused by coupling of 
the in-plane and flap blade modes at high collective pitch and is discussed 
in detail in Reference 1. 

In recognition of the importance of in-plane motion to the development 
of satisfactory hingeless rotors, the present study has investigated the 
in-plane behavior of multibladed rotors. 
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Along with the large mean hub moment capability of the hingeless rotor 
comes inevitably the capacity to produce large oscillatory hub moments. As 
long as the advance ratio stays well below p = 1.0, however, they are nor- 
mally not significant, except in certain transition flight conditions. But 
at advance ratio greater than p= 1.0, the shaft vibration force levels 
become very important. And since it appears that the next major development 
in hingeless rotor vehicles involve flying at advance ratios greater than 
P = 1.0, the present study investigated shaft-transmitted vibratory forces 
due to high advance ratio. 

A basic problem in hingeless rotor design is the creation of section 
structural properties that do not pick up excessive stresses under any 
operating or transient condition and that at the same time maintain radial 
mass and stiffness distributions which keep the desired characteristics and 
avoid instability. In the next generation of hingeless rotor aircraft, 
expected to operate at reduced rotor rpm and p> 1.0, this problem is 
especially acute. Endurance limit stresses must not be exceeded if blades 
of sufficient life expectancy are to be produced. 

Because of the fundamental nature of blade stress determination to ad- 
vanced rotor design, the azimuthal distribution of blade section flap bend- 
ing moment were investigated in this study. 

The results of these studies are expected to provide a basis for the 
systematic optimization of the design of the next generation of high speed 
compound helicopters. This includes the slowed hingeless rotor compound 
helicopter that cruises at a speed of 300 to 350 knots with the rotor 
slowed to as little as half the hovering rpm, and is characterized 
by low rotor lift and advance ratio to fx = 1.5. It also Includes, at some 
future date, the stowable rotor aircraft with the speed of fixed wing air- 
craft, low rotor lift during conversion, rotor slowed to zero rpm and then 
stowed within the fuselage. 

The specific purposes of this study were as follows: 

1 . Derive the equations of vertical and in-plane mass element motion 
for three- and four-blade rotor-body free-flight configurations with 
gyroscope-stabilized swashplates. 



2. Analyze the equations for their steady harmonically forced 
vibratory response in the form of mean and harmonic aeroelastic 
shaft-transmitted forces. 

3. Analyze the equations for blade section lift and flap bending 
moment variation with azimuth. 

4. Analyze experimental shaft vibratory, and blade load data for three 
rotors : 

a. A 33-foot 3-blade rotor with a high constant speed gyroscope- 
stabilized swashplate and fixed shaft. 

b. A 7.5-foot 4-blade rotor with fixed shaft and swashplate. 

c. The 35-foot 4-blade rotor of the XH-51A compound helicopter in 
free flight with a rotor speed gyro-stabilized swashplate and 
the shaft restrained by the free body. 

5 . Compare theory and experiment . 

6. Parameter range for which the study is valid is as follows: 


a. 

Stiff bladed rotors 

p> 

1.3. 

b. 

High advance ratio 

m > 

• 3. 

c . 

Subcritical advancing tip 

Mi 
-L r\r 

< . 85 . 

'k 


Phase I of this study (References 2, 3, and 4) investigated the vertical 
motions equations of the 33-foot 3-blade rotor. The mean hub and swashplate 
moment and thrust mean aeroelastic derivatives were calculated, neglecting 
the effects of the harmonic differential equation coefficients as is the 
practice with conventional articulated rotors at low advance ratio. Good 
agreement was obtained with experiment at low values of advance ratio. 

Cyclic pitch to trim hub and swashplate moments to zero, control effectiveness, 
and stability of the gyro-stabilized system were also investigated and 
reasonable agreement found with experiment. 

Equations of flapping motion, including aerodynamics, were found for 
the 7.5-foot rotor by Sissingh and Kuczynski in References 5 and 6. The 
mean aeroelastic derivatives of hub moment with respect to cyclic and 
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collective pitch and angle-of-attack were calculated, including effects of 
blade mean, first and second harmonic response. Good agreement was obtained 
with experiment. 

Effects of induced inflow and first and second blade flapping modes were 
investigated by Ormiston and Peters in Reference 7- This study showed the 
necessity of including more than one flap mode if the ratio of flap frequency 
to rotor rotational frequency was less than 1.3- It also showed the large 
influence of induced inflow at advance ratio M< -3 and also that significant 
effects of induced flow existed at all advance ratio. 

The necessity of including the rotor and gyroscope degrees of freedom 
in the equations of motion of the total airframe is shown by Heimbold and 
Griffith in Reference 8. 

Blade section aerodynamic lift and flap and in-plane bending moments 
presented by Bartsch in Reference 9> Volumes I and II, and comparison with 
comprehensive method results by Sweers , in Volume III, for flight tests of" 
the XH-51A compound helicopter, form a data base for assessing the validity 
of the current formulation in the modest /u , low flap frequency range. 

Deckert and Me Cloud analyzed 33-foot rotor data in Reference 10 and 
discovered that hub moments produced by cyclic pitch approach zero as rotor 
rpm and advance ratio approach zero and infinity respectively, thereby 
precluding their use in trimming mean hub moments. They also showed that 
small values of cyclic pitch were useful, at high advance ratio, in reducing 
shaft vibration. 

The experimental testing of the 33-foot 3-blade rotor analyzed under 
this contract took place in the NASA-Ames Research Center 40 by 80 foot wind 
tunnel. The maximum wind speed attained was 120 knots and advance ratios 
of M= .4, .5, .8, 1.1, and 2.0 were tested over ranges of rotor rpm. 

Four mass and stiffness configurations of the 7.5-foot 4-blade rotor 
were tested in the U. S. Army Air Mobility R and D Laboratory, Ames Direct- 
orate 7 by 10 foot wind tunnel. The advance ratio range extended from .29 
to 1.75 and flap frequency ratio extended from P = 1.22 to 2.32. 
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Flight testing of the XH- 51 A compound helicopter was performed by the 
Lockheed-California Company, Rotary Wing Division. 

The equations of motion derived in this study were programmed on the 
IBM 360 digital computer, and the CPS (Conversational Programming System) 
terminal was employed in the analysis of the experimental data. 

In the design of compound helicopters for flight to advance ratios 
greater than unity (u> 1.0) elimination of shaft-transmitted vibratory 
forces and the control of blade stress excursions are expected to be the main 
tasks. This study has indicated the nature and causes of these two phenomena 
and provided a straightforward analytic tool for their prediction. The 
groundwork has therefore been laid for the invention of devices, such as 
harmonic swashplate control, and optimization of design for the minimization 
of flight loads and vibrations in advanced hingeless rotor helicopters. 

It is also possible, with vibration and loads calculable at extreme 
values of design and flight operation parameters, that refinements in the 
methods may be made which will allow calculations in the realm of the con- 
ventional helicopter. 
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SYMBOLS 


Symbols and matrices employed in the vertical motions equations are 
shown first. Those required for the in-plane equations follow, and sub- 
scripts common to both systems complete the section. 


b 

b .m. 

B 

c 

c i 

C 

m 


0 


e 


c 

m 


e 


c 


l 


4 > 


c 


T 


C t 


a 


C R 


c 


s 


D 


F, t 


number of blades 

blade flap bending moment ft -lb 

tip loss factor 

blade chord ft 

blade root cutout fraction 


hub pitch moment coefficient 

hub roll moment coefficient 
s-washplate pitch moment coefficient 

swashplate roll moment coefficient 

thrust coefficient 

blade section lift curve slope 
feathering friction 
swashplate damping 
drag 

jet engine thrust 


M 

p(UR) 2 ttIO 

L 

p OR) 2 ttR^ 
M 6 

p(uR) ? rrR 3 

M 4 > 

pOr) ? ttr3 


p(QR) ? ttR 2 


ft -lb /rad/ sec 

ft-lb/rad/ sec 

lb 

lb 
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blade flap generalized force 
rotor flap generalized force 


blade flap moment of inertia 

blade moment of inertia about 
the quarter chord 

gyroscope diametral moment of 
inertia 

airframe pitch and roll moments 
of inertia 

mechanical advantage blades to 
swashplate 

swashplate spring 
lift 

hub rolling moment 
hub pitch moment 
blade mass distribution 
rotor flap generalized mass 


swashplate pitch moment 


swashplate roll moment 

n^ flapping mode 
th 

p blade 

ratio of blade first flap frequency 
to rotor rotation rate 

blade pitch moment about the quarter 
chord 


2 

slugs ft 

2 

slugs ft 

2 

slugs ft 

2 

slugs ft 

ft -lb /rad 
lb 

ft -lb 
ft -lb 
slugs/ft 

ft -lb 
ft -lb 


ft -lb 



q 

dynamic pressure 

lb/ft 2 


r 

radial distance from shaft centerline 

ft 


R 

rotor radius 

ft 


t 

time 

seconds 


T 

thrust 

lb 


V 

forward speed 

ft/ sec 


v 

blade root shear 

lb 


w 

airframe gross weight 

lb 


x, y> z 

rotor rotating axes 



x, Y, z 

stationary axes 



a 

angle-of-attack 

deg or 

rad 

P 

blade linear flap angle 

deg or 

rad 

P. 

o 

rotor precone 

deg or 

rad 

y 

fraction of critical damping 

4 

Cq cR 

p ot 

y 

blade Lock number 

h 


h 

blade parabolic mode tip deflection 

ft 


b o ,b e ’ 6 d 

e 

rotor flapping mode deflection 
swashplate tilt 

ft 

rad 


e o 

collective pitch 

deg or 

rad 

0 .75R 

collective pitch at three quarter 
radius 

deg or 

rad 

e 

t 

blade twist rate 

rad/ ft 


®c/U (or 9p) 

blade pitch about the quarter chord 

rad 


6r 

blade feathering pitch 

rad 




deg or rad 


e 


lc, 


'Is 


® 

A 

F 

P 

cr 

<t> 

<t> 

4^ 

w 6 


SI 


SI 


G 


M 

[ a r] 

[af] 

[B] 

[ b r] 

[cf] 


cyclic pitch 

rot or -air frame pitch angle 
b],ade forward sweep 
advance ratio 
air density 
solidity 

swashplate roll tilt 

rotor-airframe roll angle 

azimuth position of rotor 

position of p th blade relative to 
blade 1 

cant angle 

blade first flap frequency 

rotor rotation rate 
gyroscope rotation rate 


rad 

deg rad 


slugs/ ft 


3 


rad 

rad 

deg or rad 


rad 


deg or rad 
rad/ sec 

rad/ sec 
rad/ sec 


Matrices 


aerodynamic stiffness, rotating axes 
aerodynamic damping, rotating axes 

aerodynamic forcing, rotating axes 
aerodynamic stiffness, stationary axes 
aerodynamic damping, stationary axes 
centrifugal and structural 
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Cd 

[Da] 

DO 

P] 




M 

IP} 

H 


transformation, rotor to blade freedoms 
damping 

centrifugal and structural, stationary axes 
forcing matrix, stationary axes 

rate of change of blade force with blade 
deflection 

inertia 

airframe inertia 

airframe aerodynamic stiffness 

airframe aerodynamic damping 

airframe acceleration aerodynamics 

centrifugal and structural, stationary axes 

transformation from stationary to rotating 
axes 

forcing angles 

rotor degrees of freedom 

blade degrees of freedom 


Additional symbols employed in the in-plane motions equations 


a.f- blade axial force lb 

c ^ minimum drag coefficient 

o 


e 


pivot radial position 


ft 



F ip 


h • H 

o u 




*b 


shaft 


Divot 


R 

V V 


; , k 

e b % 


“r 


n. f . 


N 

p • t • 

r 

c.g. 


s . t . 


u 

V 

X 

X 


blade in-plane force lb 

rotor in-plane mode generalized force 

/■ 

blade moment of inertia about the shaft slugs ft 

blade moment of inertia about the pivot slugs ft 

hub moment of inertia 
blade inertia parameters 


slugs ft 


body lateral spring to ground lb/ft 

body longitudinal spring to ground Ib/ft 

lift on p^ blade lb 

blade mass slugs 

hub mass slugs 

blade normal force lb 

rotor yawing moment ft -lb 

blade pivot torsion ft-lb 

blade center of gravity radial ft 

distance 

blade shaft torsion “ ft-lb 

blade normal displacement ft 

blade axial displacement ft 

longitudinal displacement ft 

longitudinal force lb 
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y 

Y 


C 

&o' W ^n’ Cd 




lateral displacement 
lateral force 

blade rotation about its pivot 
rotor mode deflection 

ratio of inclination of blade 
in-plane principle axis to 
feathering pitch 

blade rotation about the shaft 


uip 


[cor] 


blade non-rotating undamped natural 
frequency 

Coriolis force matrix 


Subscripts 


0, lc, Is, 2c, 2s . . . . harmonic component 

c/4 about the quarter -chord 

d differential collective mode 

f feathering 

g gust 

i lift 


n 

n 

P 

1 , ? 
r 


normal to blade leading edge 

n^* 1 mode or harmonic 
th 

p blade 

p , ♦ . . , b blade numbering sequence 
rotating axis system 


ft 

lb 

rad 

rad 

rad 

rad/ 


• 75P 


blade three-quarter radius 



xy 

XY 

z 


in-plane equations op in-plane 
rotating axes 

in-plane stationary axis system 

vertical equations or vertical system 
rotating axes 

vertical system, stationary axes 


Less frequently used symbols are defined in the text. 



HINGELESS ROTOR TYPES STUDIED 


In this report a set of stiff -bladed hingeless rotors operating at high 
advance ratio is considered. Its selection was based on the desire to study 
rotors that could be employed in future vehicle concepts. Therefore, those 
aspects of the advanced helicopter design that cross into the well-understood 
areas of conventional helicopters, such as hover and low advance ratio flight, 
are not studied. 

This philosophy crystallized Into the consideration of rotors which could 
be applied to the two following types of vehicle: 

1. Slowed hingeless rotor winged compound helicopter 

2 . Stowable rotor aircraft 

These considerations resulted in the decision to study only the high advance 
ratio flight range and rotor blades with flapping stiffness in the 

range 1.3<P< oa * It further resulted in limiting the flight regions to low 
lift values and subcritical advancing blade tip Mach numbers. 

The blades, experimentally investigated, had Inoptimum structure and there- 
fore low Lock numbers, but this is not considered to restrict the applicability 
of the theory correlation. 

Because the rotors were of the hingeless type, the minimum number of 
blades considered was three. Four-blade rotors were also extensively studied. 
The theoretical derivations were extended to any number of blades. 

Since such rotors can be employed with gyroscope- stabilized feathering 

systems, blades swept forward of the feathering hinge and corresponding feath- 
ering moments and degrees of freedom were considered throughout the analyses. 

Geometry Variations 

The variations in blade geometry considered in this study are as follows: 

R ft 

c ft 


Radius 

Chord 

Sweep angle 


A 


radians 
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Cutout fraction 


Precone 
Twist rate 


radians 

radians/ft 


These are illustrated in Figure 1. It should be noted that both the 
locus of quarter chord points (along the r-axis) and the feathering axis 
pass through the center of the shaft or center of rotation (the z-axis). 

In blade configurations in which the quarter chord or hinge lines do not pass 
exactly through the center of rotation, fictitious hinge lines and quarter 
chords may be assigned with little loss of precision. 

The analyses presented in this report are valid for precone and twist 
values that do not seriously violate small perturbation values (say up to 10° ) 
and for any value of blade radius. Chord lengths, on the other hand, should 
not exceed a chord-to-radius ratio of about l/lO or blade section aerodynamic 
pitch damping, neglected in the theory, may become significant. Blade tor- 
sional deformation, neglected in these analyses, can also result in significant 
changes in blade airloads should the ratio of aerodynamic moment to torsional 
rigidity become large enough. 

Sweep angle and root cutout values may be of any magnitude . 



Figure 1. Blade Geometry 




The particular equations derived to completion in this report are for 
three- and four-blade rotors but indications are given at every step as to 
the procedure to follow to extend the derivation to rotors with 5, 6,.... 
n blades. 

In considering the complete rotor, an axis system rotating with the 
blades is required to act as a reference for the location of each of the 
blades. It is denoted by lower case letters. The axis system is selected 
so that the positive z-axis runs along the shaft upward, in its nominal un- 
tilted position. The axes rotate with the rotor but are not attached to it, 
for the rotor can pitch, roll, plunge, surge, sideslip, and yaw relative to 
the rotating axis system. The x-axis runs along the nominal position of the 
number one blade quarter-chord locus, if precone is zero (otherwise normal to 
the z-axis). The y-axis extends normal to the x, z plane, positive to the 
right, to form a right-handed axis system. 

The blades are numbered in ascending order moving around the rotor in a 
clockwise direction so that, to an observer on the ground watching the blades 
go by, the blades would appear in the order 1, 2, n, 1, 2, 3*** 5 

The arrangement of blades and rotating (lowercase) axes is shown in Figure 2. 



Figure 2. Rotating Axis System and Blade Numbering 


17 





The figure also shows schematically the attachment of the blade feathering 
horns to the swashplate via the stabilization gyroscope, in the general case of 
this report. Tilting the swashplate, and the parallel gyroscope, causes the 
blades to feather cyclically. The relationship between swashplate tilt and 
blade feathering is shown in Figure 3. 

Swashplate tilt relative to rotating axes is denoted by 0 , 0^ pitch 

and roll about the y and x axes respectively. Because of the blade feath- 
ering arm and linkage geometry, pitch and roll do not separately produce pure 
"cyclic pitch” of the blades as is usually the case with articulated rotors, 
("cyclic pitch" here is relative to rotating axes.) This slewing around of 
the gyroscope (and swashplate) relative to the blades is denoted by i|f , the 
cant angle, and is shown in Figure 3- By this definition, most articulated 
rotors would possess a cant angle of ninety degrees (\|t o = 90 °). 
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Cant angles different from 90° are employed in gyroscope -stabilized 
hingeless rotors for purposes of free-flight stability. If the shaft is 
prevented from pitching and rolling, as is the case during wind tunnel 
testing, the rotor-gyroscope stability is not affected by cant angle. It 
should be noted that the cant angle is relative to the blade feathering axis. 
The blade quarter chord leads the feathering axis by the sweep angle. 

In order to allow the use of a smaller gyroscope the gyro arm is made 
shorter than the blade feathering horn arm so that the gyro tilt angle is 
somewhat greater than the maximum feathering displacement. This is called 
mechanical advantage . 

. Horn arm n A , 

k = = Mechanical Advantage 

Gyro firm 

Mechanical advantage is such that values greater than about 1.5 are 
difficult to achieve mechanically. 


Degrees of Freedom 

The geometric form of the rotor systems studied in this report has been 
discussed in the previous section. It was necessary to establish a reference 
axis system rotating with the rotor in order to describe the blades, linkages, 
gyroscope, and rotating part of the swashplate. Similarly, in order to dis- 
cuss the rigid body to which the rotor is attached and the stationary part of 
the swashplate, it is necessary to establish a second reference axis system 
that is stationary. Denoted by capital letters, this system employs the same 
vertical or Z axis as the rotating system. Its X and Y axes, however, 
do not rotate but remain fixed relative to the earth. (in actual fact they 
translate with the mean forward speed of the vehicle.) The X or longitudinal 
axis is positive in the nominal aft direction of the body and the Y axis ex- 
tends normal to the X, Z plane positive to the right, again forming a right 
handed system. Both axis systems are shown in Figure 4. 

The rotating axes are related to the stationary axes by the rotation 
angle ty, which is a function of time. Since the rotation rate is assumed 
to be constant in these analyses, i|r = fit 
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Figure 4 . Stationary and Rotating Axes 


where fl = the rotor rotation rate - radians/sec 
t = time in seconds. 

Perturbational motions of the rotor and shaft, described in terms of 

degrees of freedom in the rotating axis system, may be transformed to station- 
ary axes. In this form the shaft motions are in a convenient form to Include 
the rigid motions of the vehicle body. 

It is clear, therefore, that the degrees of freedom of the system may 
be described in rotating or stationary axes. (There are some exceptions which 
will be noted later.) Motions of rotating parts are basically described rela- 
tive to rotating axes and motions of stationary parts relative to stationary 
axes. The complete set of motions equations could be described in either 
coordinate system. In these analyses, however, the rotor alone is described 
in rotating axes; its equations of motion are transformed to stationary axes 
and these are combined with the body equations of motion and solved in 
stationary axes. 
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Whether thought of relative to rotating or stationary coordinates, the 
rotor-gyroscope-body motions have been separated ;yixo two sets for the pur- 
pose of this study; those that produce perturbational mot^§is of the elements 
of mass essentially in the vertical direction, and those that produce mass 
element motions essentially in the plane of the rotor disk. The two sets of 
degrees of freedom can only be separated in the small perturbation motion 
sense since nonlinear coupling between the two sets becomes significant with 
increasing blade loads and displacements. 

It should be noted that the vertical motions equations are solved com- 
pletely independently of the in-plane motions equations. The in-plane motion 
equations forcing functions, however, are produced almost entirely as a by- 
product of the motions of the vertical equations degrees of freedom and their 
aerodynamic forcing functions. 

The two sets of degrees of freedom will be described separately in the 
following: 

Vertical motions of mass elements * - The degrees of freedom employed in 
describing those motions of the rotor-gyroscope-body system that are essen- 
tially in the vertical direction can be most easily pictured relative to sta- 
tionary axes (with one exception, discussed later). Figure 5 shows the sta- 
tionary X, Y, Z axes degrees of freedom of a 5-blade rotor. The set for 
3-blade rotors is obtained by deleting the two scalloped disk flapping modes 
6^c and 6^ s , the reactionless flapping cervine and sine degrees of freedom. 

For rotors with even numbers of blades 4, 6, etc, in addition to the 

reactionless modes existing in pairs, there is a single reactionless mode, 
sometimes called the differential collective mode (or degree of freedom). It 
is characterized by the blades being alternately up and down. 
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REACTIONLESS FLAPPING 
COSINE COMP. 



REACTIONLESS FLAPPING 
SINE COMP. 




In this report the degrees of freedom are arranged in vectors in a fixed 
order to facilitate matrix algebra. For 3- ? 4-, and 5-blade rotors the 
vectors of degrees of freedom in stationary axes are as follows: 


3-blade rotor 4-blade rotor 5-blade rotor 



NOTE: 1. Subscript Z denotes stationary axes. 

2. Subscript d denotes differential collective mode. 

The degrees of freedom have been shown in stationary coordinates for 
the sake of clarity but they, in fact, represent the end product; the form 
in which they are employed in the equations of motion to be solved. 

In order to obtain these degrees of freedom a parallel set in rotating 
coordinates must be transformed to stationary axes. The rotating axis de- 
grees of freedom employ the same symbols as the stationary axis modes but are 
differentiated by subscript "r." They look the same as the stationary modes 
but the mode shapes are relative to the rotating x, y, z axis system and 
remain undistorted as the rotor rotates or are independent of azimuth (the 
number one blade lies along the x-axis at all values of azimuthal orienta- 
tion of the rotor). 

The transformation from rotating axis degrees of freedom to stationary 
axis degrees of freedom must therefore be a function of azimuth or time if 
the rotation rate is constant. 
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The transformation for the 5-blade rotor is as follows: 



z - rotating axes 
Z - stationary axes 
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For the 3-hlade rotor the transformation is as above with the last two 
rows and columns deleted. For the 4-blade rotor the last row and column are 
deleted and the lower right corner element cos 2 t|r is replaced by 1. 

For aircraft body and gyroscope (or swashplate) angular motion and rotor 
disk tilt, unit values of displacement in stationary coordinates are equivalent 
to displacements in rotating coordinates which vary sinusoidally once per rev- 
olution with azimuth. For example, a unit value of body pitch in stationary 
coordinates 0 = 1.0 radian is equivalent to the following sinusoidal varia- 
tions with azimuth in rotating coordinates: 

© r = 0 cos t 

$ = 0sin f 

r 

Body plunge, collective flapping and, for the four-blade rotor, the 
differential collective mode are the same in both coordinate systems. 
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Constant reactionless flapping displacement in stationary coordinates, 
however, is equivalent to reactionless flapping in rotating coordinates that 
varies with the second harmonic of azimuth. For example, if the disk of a 5- 
blade rotor took up a constant cosine scalloped shape in stationary coordinates, 
f> 2 c ,the equivalent motion in rotating coordinates varies with azimuth as follows 


6 


2c 


r 


6 


2s 


r 


& 2 c cos 2 t|i 
§ 2 c sin 2 i|r 


The relationship between rotating and stationary axes vectors of degrees 
of freedom is also expressed in this report for the sake of conciseness, in 
vector notation, as follows: 



z = rotating axes 
Z = stationary axes 


Although the equations are employed in stationary coordinates, they came 
into that form through the application of the - j T z J transformation to the 
more fundamental equations of motion derived relative to rotating coordinates. 
It is necessary, therefore, to first develop the equations relative to ro- 
tating coordinates or the relationship between the motions and the generalized 
forces of the rotating degrees of freedom. 

The first step in this development is the determination of the kinematic 
relationship between the motions of the individual blades in their single- 
blade degrees of freedom and the motions of the complete rotor in its rotating 
degrees of freedom. The relationship is independent of the azimuthal position 
of the rotor. This forms the subject of a later section. 

Motions of the rotor in its rotating degrees of freedom may be described 
in terms of the single blade degrees of freedom (shown in Figure 6). 

The vertical single blade degrees of freedom may be arranged in a vector 
to facilitate matrix algebra. 
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VERTICAL MOTION 



THEREFORE RELATIVE FLOW IS DOWNWASH 



PARABOLIC FLAPPING MOTION 


<b ROTATION 



6,5 TIP UP + 


Figure 6. Single Blade Vertical Degrees of Freedom 


29 




N 

The pitch, roll, and plunge motions of the rotor disk are composed of the 
first three blade freedoms: linear flap P, vertical plunge z, and blade 

feathering 9 p . Gyroscope (and swashplate) pitch and roll are composed of 
blade feathering and linear flapping (due to the sweep of the feathering hinge). 
Rotor disk collective, pitch and roll and reactionless flapping modes are all 
defined in terms of the blade parabolic or first flap structural mode. 

The vertical motions equations contain only the first blade dynamic 
flapping mode . 

In-plane motions of mass elements . - The stationary axis degrees of free- 
dom employed to describe rotor motions essentially in the plane of the disk 
are as shown in Figure 7 . The freedoms shown are for a five-blade rotor. 
Deletion of the cosine and sine reactionless lead-lag freedoms provides the 
set needed for a three-blade rotor. Four -blade rotors require the addition 
of one differential collective lead- lag mode to the three-blade set. It is 
the same in rotating and stationary axes. 

The number of blades shown in Figure 7 is arbitrary; the displacements 
merely show the deflected positions the blades would occupy at that azimuth. 

It should be noted that positive lead-lag is lead or counterclockwise structural 
deflections. (This is opposite the convention in Reference 11 but is used 
since all counterclockwise displacements are to be positive, as is the normal 
mathematical convention.) At this point, the blade in-plane deflection mode 
has not been defined. 

The' in-plane degrees of freedom in stationary axes are arranged in vectors 
to facilitate matrix algebra, as were the vertical motions freedoms. 
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Figure 7. In-Plane Motions Degrees of Freedom in Stationary Axes 
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TORSION LEAD-LAG ? 
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Figure 7. (Continued) 
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Figure 7. (Concluded) 
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3-blade rotor 4-blade rotor 5-blade rotor 



Compared to the vertical motions degrees of freedom, a motion equivalent 
to swashplate tilt is missing* The in-plane motions are therefore described 
by two fewer degrees of freedom than the vertical motions. 

The next step is to visualize these modes in axes rotating with the 
rotor. The lower case x, y label for rotating axes replaces the uppercase 
X, Y of the stationary axes system but the modes maintain the same shape. 

The major difference now is that the blades maintain fixed azimuth positions 
relative to the x, y axes. The x-axis is always directed through the axis 
of rotation parallel to the number (l) blade quarter chord line and the blades 
are numbered in ascending order as they pass a fixed azimuth position. The 
azimuth location of the pth blade relative to number (l) blade is 

♦ P = - (p - 1} m 

where b is the number of blades. 

The rotating axis degrees of freedom employ the same symbols as the sta- 
tionary axis freedom but differentiated by subscript "r. fl 

Vectors describing the in-plane rotor motions in rotating degrees of 
freedom may be transformed to stationary axes degrees of freedom by a sine- 
cosine transformation very similar to the one employed for the vertical mo- 
tions. For five-blade rotors, the transformation of rotating axes freedoms 
to stationary axes is as follows: 
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For 3-blade rotors, the last two rows and columns are deleted. For 
4-blade rotors, 1.0 is added to the lower right corner of the 3-blade mo- 
tions expanded by one row and column. 

The above relationship is also expressed in matrix form as follows: 


xy 


= T 


xy| XY 


Vector subscripts 
xy - rotating axes 
XY - stationary axes 


Before the equations of motion in rotating axes can be developed it is 
necessary to determine the kinematic relationship between motions of indi- 
vidual blades in their single-blade degrees of freedom and the motions of the 
complete rotor in its degrees of freedom relative to axes rotating with the 
rotor. This is the subject of a later section. 


Motions of the rotor in its rotating degrees of freedom may be described 
in terms of the single blade degrees of freedom shown in Figure 8. The first 
flexible mode has been defined to be a straight line lead-lag with an offset 
pivot. This mode shape was assumed after an inspection of experimental mode 
shapes for the 33-foot rotor. Examining the reasons for its very close ap- 
proximation to this shape, it became evident that all hingeless rotors with 
relatively flexible root sections in-plane and relatively stiff in-plane blade 
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RIGID BODY LEAD-LAG TRANSLATION u 



RIGID BODY SPANWISE TRANSLATION v 





Figure 8. Single Blade In-Plane Degrees of Freedom 





stiffness outboard would have shapes of this type. The specific shapes and 
frequencies , etc., will be discussed in a later section, as they vary with 
rpm. 

The single blade degrees of freedom may be arranged into a vector for 
matrix operations: 



The rotor in-plane degrees of freedom relative to coordinates rotating 
with the rotor are composed of the single blade degrees of freedom. The only 
structural dynamic mode included is the first in-plane lead-lag mode denoted 
£. Higher dynamic modes have not been included in this analysis. 

The lateral and longitudinal degrees of freedom y^ and in ro- 

tating coordinates are composed only of u and v single blade degrees of 
freedom. The rotational mode is composed only of the blade rigid body 

rotation mode y. All the remaining rotor degrees of freedom, relative to 
rotating axes are lead-lag modes and are composed solely of flexible first 
lead-lag blade degrees of freedom £. 

Kinematic Relationships 

In this section the rotor degrees of freedom relative to rotating axes 
are defined in terms of single blade degrees of freedom. The kinematic re- 
lationships between the two sets of freedoms, for each set of motions, verti- 
cal and in-plane, are expressed as transformation matrices 

Sets of degrees of freedom for vertical and in-plane motions are treated 
separately. 

Vertical motions blade-rotor kinematics . - Arranging the single blade 
degrees of freedom in a vector facilitates the matrix algebra and allows the 
displacements of the pth blade to be represented by a single symbol T| 

z 
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where 



6 

P 


The blade motions may then be expressed conveniently in terms of rotor 
motions in degrees of freedom relative to rotating axes. 



where p indicates the blade number. Rewriting: 
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for a five -blade rotor 
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The 

next task is 

to define 

the 

motions of the 

individual blades in terms 

of those 

of the 

rotor 

, the elements 

of the matrix. 

The motions of 

the blades 

are relative to 

the blade line 

of aerodynamic 

centers (or 

quarter 

chord) 


projected Into a plane normal to the undeflected rotor shaft. 


Consider a rotor with b number of blades with the number 1 blade located 
along the fx axis (the axes rotate with the rotor). The position of the 
locus of aerodynamic centers of blade p is located at: 



x 
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Blade motion due to shaft motion - Blade rigid flapping and pitching about 
the aerodynamic center are the only blade motions resulting from shaft pitch 
and roll displacements. 



I 

x 


The rigid flapping of an arbitrary blade p due to shaft motion is: 

B = - © cos ‘if + $ sin i|r 
p r T p r T p 

/ \ 2 TT , . / \ 2 TT 

9 = - ® cos (p-1) -r— - * sm (p-1) ”T— 

Kp r ^ r b 

The pitching of blade p due to shaft motion is: 

6 = ® sin \lr + $ cos ilr 

p r p r p 

- . /nN 2n , f x 2 tt 

9 = - ® sm (p-1) r + * cos (p-1) — r— 

p r v ^ ' b r v ' b 

Rotor plunge motion produces only blade plunge motion. For a blade p ? 
the relationship is merely: = z. 


In matrix notation the blade motion due to shaft motion is: 
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where p = 1, . . , b. 


Blade motion due to gyro motion: Gyro pitching and rolling motion causes 

blade pitching and flapping relative to the blade quarter chord axis. In 
order to resolve the gyro motions into blade motions, define the axes x, and 
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Positive values lor 6^ and 0 b swashplate deflections produce 
positive blade feathering angles at i|) = l80* and 90° respectively. Note 
that the blade feathering axis is located behind the quarter chord axis by 
the angle A. However, the azimuth position is measured relative to the 
quarter chord axis. 

With the above in mind, the feathering angle of a blade located at 

is : 

0 f = 5 \ cos * p + \ sin *p] 

e f - - 1 1 c ° s (p - 15 y - 1 ^ sin ‘p- 1 * y 

P 


or in matrix notation: 

\\\‘ ' 5 t° S <p ' 1) y 



where k is the gyro to blades mechanical advantage. 


The relationship between blade feathering about the feathering axis to 
pitching and flapping relative to the quarter chord axis is: 
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Now the motions about the quarter chord can be expressed in terms of the 


swashplate angles 0^ and 0^: 
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therefore : 
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Blade motion due to rotor flexible flapping motion: Blade flapping 

displacements are limited to the first flap mode, various combinations of 
which make up the rotor flapping modes. In rotating coordinates, the number 
of rotor modes will then equal the number of blades. A 2-blade rotor would 
have a collective and a pitch (or roll) flapping mode; whereas a 3-blade 
rotor would have both pitch and roll flapping modes, plus the collective 
mode. A 4-blade rotor would additionally have a self-balancing rotor mode, 
which would produce no inertial forces in the earth-fixed axes system. A 
5-blade rotor would have two self-balancing modes. 


For example, the modes of a 4-blade rotor would be as follows: 
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The relationship between the rotor flapping modes and the flapping of a 


blade at ib is: 
Y V 


6_ = - 6 nQ cos n i|;„ + sin n t 




6 = - 6 0 
p n6 


cos n (p-l) - 6 . sin n (p-l) 

r n ^r 


For rotors with an even number of blades: 

n — 1 ^ • • • • j 0 • 5b 

For rotors with an odd number of blades: 
n = 1,...., 0.5 (b-l) . 
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where n is a number representing the harmonic shape of the mode relative 
to the axes rotating with the rotor. For example, n = 1 represents rotor 
disk tilting relative to the shaft; n = 2 represents a scalloped shape for 
the rotor tips with two lobes up and two down regardless of the number of 
blades . 

Substituting the required values of n yields the modes. Thus: 

6 p = ■ 6 e cos ( p_1 ) ur " 6 d sin ( p-1 ) or n = i 

y r 

6 = - 6 cos 2 (p-l) 1— - 6 sin 2 (p-l) 1— n = 2 

p pc b 2s ^ b 
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In addition, the collective mode will be: 

6=6 
P o 


In matrix form: 



In summary, the transformation matrix relating blade motions of the 
degrees of freedom in rotating axes are as shown on the following page for the 
five-blade rotor: 


In-plane motions blade-rotor kinematics - With the in-plane displace- 
ments of the pth blade represented by the vector: 


^5 
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the relationship between the dispacements of degrees of freedom of the rotor 
relative to rotating axes and the displacements of the pth blade may be ex- 
pressed as follows: 




where B is the vector of displacements of the rotor degrees of freedom 


relative to rotating axes. 


Determination of the elements of the f D j matrix, or the pth blade mo- 

[ xy pJ 

tions that occur due to the motions of each of the rotor degrees of freedom 

follows . 
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DUE TO UNIT x DISPLACEMENT: 


FOR THE P TH BLADE OF A 
"b" BLADED ROTOR 



48 



DUE TO UNIT f DISPLACEMENT: 

o 




1.0 


e to unit Q displacement: 

C = - cos \|r 
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e to unit C x displacement: 

C = sin \|r 
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e to unit displacement : 
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Due to unit C 2s displacement: 
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If more than five blades are employed in the rotor design, then 
eral expression for blade displacement due to unit reactionless mode 
ments can be used. 

For unit reactionless mode dispacements £ and Q , where m 
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2 . . .n 
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and n = for odd numbers of blades 
b = number of blades 

the displacements of the pth blade become 

C = - cos mil; 

P P 

C = sin mil; 

P P 

Theoretical Considerations 

Some of the physical elements of the rotor-gyroscope-body system or the 
forces on, or motions of them, were simplified for the sake of clarifying the 
behavior of the complete system. The chosen simplifications did not greatly 
distort the completed result and the reasons for this are discussed in this 
section. 

Those aspects of the system elements simplified included: rotor aero- 

dynamics, blade section aerodynamics, radial flow effects at blade tips, in- 
plane components of aerodynamic lift forces, in-plane motion of a single blade, 
blade flap to in-plane inertia coupling, and the high-speed gyroscope physical 
representation. The ranges of system parameters over which the simplifications 
are expected to be valid are indicated. 

Some concepts useful in understanding the behavior or hingeless rotors 
are also discussed. These include mean rotor aeroelastic derivatives, the 
composition of the hub force two-component vibration derivatives, residual 
forces and how these may be trimmed by the application of cyclic pitch. 

Aerodynamic Simplifications . - Flight at advance ratio greater than 
|jl = 0.7 has not been studied very extensively, up to this time. Therefore, 
the present study should be considered exploratory in nature. The simplest 
concepts which yield approximately correct answers have been used. 

In this study of the basic behavior of hingeless rotor systems, the math- 
ematical models representing them have been shorn of all but the most essential 
considerations. Some of the effects eliminated would have contributed signifi- 
cantly to the magnitudes of the final answers at certain critical conditions, 
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but including them would have obscured the fundamental behavior of the system 
(and more comprehensive methods are available for the purpose of producing 
accurate quantitative results). 

Some of the simplifications are as follows: 

Linear aerodynamics: The present tests have been limited to conditions 

which do not seriously violate aerodynamic linearity. Allowing the use of 
linear aerodynamic theory opens the door to straightforward techniques for 
estimating such things as radial lift induction effects, effects of downwash 
from previous blade passages and unsteady aerodynamics. Only one of these 
effects has been investigated in this study. The main consequence, however, 
of aerodynamic linearity is that superposition of the effects of angle changes, 
for example: control angle, angle of attack, and precone may be calculated 

* independently of each other and summed up to yield the correct answer. 

Aerodynamics is expected, to be linear as lour as tne combined sec lion 
angle-of-attack remains less than, say, 10 degrees. Luring the conversion 
phase of flight of a stowable rotor aircraft and high-speed flight of slowed 
rotor compound helicopters, the nominal rotor lift is small so that blade 
angles -of- attack are well below stall limits in areas of significant dynamic 
pressure. In addition, over most of the conversion, tip speeds are well below 
the speed of sound. 

In the formulation of the equations of motion the inertia and structural 
forces are also kept within the linear range,- except for in-plane forcing 
functions, so that the theory of systems of linear ordinary differential 
equations, with harmonic coefficients, can be employed. 

Rotor-induced inflow: At high advance ratio, if induced inflow is 

assumed uniform over the disk, its value is very small. Even though uniform 
inflow is generally considered a poor approximation, it is expected that at 
very high advance ratio even the nonuniform induced inflow effects are neg- 
ligible;- therefore, it has been completely ignored in this report. 

Advance ratio effects: With any forward speed at all there is a region 

just to the left of the rotor mast where the net flow proceeds past the air- 
foil sections from the trailing edge to the leading edge. It is small at 


52 



low advance ratio and approaches 50 % of the disk as the advance ratio 
approaches infinity. Effects of reverse velocity have been explicitly ac- 
counted for. The blade trailing edge has become the new "leading edge" and 
section aerodynamic centers have been assumed shifted to the three quarter 
chord locus and local lift curve slope has been reduced. 

Another ramification of advance ratio is that it controls the geometry 
of the vortex patterns shed by the individual blades. At low advance ratio 
the vortex structure and attendant downwash from many blade passages accumu- 
lates over the rotor disk and drifts slowly downstream causing large in- 
duced downwash over the disk concentrating toward its aft edge. However, at 
high advance ratio, the tip vortices stream almost straight downstream from 
the blade tips so that downwash cannot accumulate. This is why induced down- 
wash at high advance ratio has been assumed to be negligible. 

At azimuths remote from 90 and 270 degrees at high advance ratio, the 
flow approaches the blade obliquely. In the disk aft semicircle the flow 
meets the blade flowing obliquely outboard. In the disk forward semicircle, 
it meets the blade flowing obliquely inboard. Wind tunnel tests have shown 
that resolving the flow into components parallel to the blade and normal to 
It, ignoring the radial component, and treating the blade section immersed in 
the normal component as two-dimensional flow, yields a close approximation to 
the measured pressure distribution and lift at the section (Reference 12). 
This principle has been used extensively for reducing the drag of wings of 
airplanes flying at transonic Mach numbers. In addition, recent tests have 
shown that the primary effect of spanwise flow is to increase the maximum 
lift before stall on the section and not to change the lift curve slope or 
linearity (Reference 13). 

At high advance ratio, radial flow acting on the tips of slender blades 
at a local angle-of~ attack, creates some lift. If the rotor angle-of- attack 
were zero, the blades in the forward semi-disk would still present a local 
angle-of-attack due to precone. A simplified analysis showed that even at an 
advance ratio as high as jjl = 2.0 the effects on aeroelastic derivatives of 
blade tip lift were dwarfed by blade loads induced by flow normal to the 
leading edge. Blade tip slender-body lift due to radial flow was therefore 
not included in these analyses. 
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Unsteady aerodynamics: The primary effect of unsteady aerodynamics can 

be seen by examining the growth of lift on a blade section after a sudden 
change in angle- of- attack . Lift growth is given by Wagner in terms of numbers 
of chords travelled (Reference It). Lift starts at 50$ of maximum and grows to 
90$ within 6 chords travelled. In a rotor with blades of small chord the lift 
becomes 90$ of steady state in a small fraction of the rotor tip perimeter re- 
sulting in effectively steady-state conditions being reached within less than 
10$ of total azimuthal travel if the forward speed were zero. At high advance 
ratio the tip velocity at ^ = 90° is greater than QR so the lift would grow 
to steady state in somewhat smaller azimuthal displacement. In the reverse 
flow region, the opposite would be true and a much longer azimuthal travel 
would be required to attain a near steady-state condition. Unsteady effects 
would cause some deviation from the results expected with steady aerodynamics, 
especially in the reverse flow region but there the low dynamic pressure makes 
the lift sensitivity small so the effect would be lost in the total aerodynamic 
derivative. The net effect of unsteady aerodynamics is not expected to reduce 
the effective blade section lift-curve slope by more than a few percent. For 
this reason unsteady aerodynamics is not expected to change the basic charac- 
teristics of solutions found without it and was not included in these analyses. 

Blade lift radial induction effects: - The most significant effect of 
radial induction is loss of lift near the blade tip. This effect depends on 
the aspect ratio of the lifting surface. The blades of typical rotors have 
aspect ratios greater than ten. For such high aspect ratios the tip effect is 
hardly noticeable inboard of about 80$ of the radius, and can be approximately 
accounted for by assuming the blade to be smaller in radius by a factor B 
(B « 0.97 in these analyses) than it is. 

A second effect is a small reduction in lift over the whole blade. This 
may be accounted for by reducing the local lift curve slope slightly. A re- 
duction factor of .95 was found to be appropriate for the 33-ft rotor. The 
lift at a section with the above corrections depends only on the local dynamic 
pressure normal to the leading edge and the normal angle-of-attack at the sec- 
tion. This is also known as strip theory. 


5 ^ 



In-plane component of blade lift: In the formulation of the in-plane 

equations of motion of the blades-hub mass system, it is assumed that all mass 
motions are in the plane of the disk. In actual fact the principal axis of 
the blade in-plane motions of the three-quarter radius section, for ex- 
ample, are not normal to and parallel to the disk plane but are at a small 
angle to the reference plane . 

It is, therefore, important to calculate the component of section aero- 
dynamic lift in the direction of the in-plane principal axis even though it 
continuously changes direction with azimuth. So long as these in-plane equa- 
tion external forces are calculated in the actual direction of the in-plane 
modes they may be applied to the planar equations without introducing 
significant error. 

The blade section lift force is predominantly in the vertical direction 
so that the in-plane component magnitude is controlled by the sine of the 
small angle between the in-plane principal axis and the relative wind vecotr 
component normal to the blade leading edge. This assumes that the blade sec- 
tion resultant force is perpendicular to the wind velocity vecotr and may be 
justified by examining a section of a blade with local angle of attack applied 
to it. The lift and drag are available from section theory and experiment 
(Reference 12). 




The drag is very small and parallel to the wind vector, it is at c^ = 1.0, 
for example, 2 % of the lift or c^ = .02. This causes about a 1.0* change in 
the resultant vector from the position of the lift vector alone. If the 

o 
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or drag at zero lift is considered separately, then the resultant is less 
than l/2 degree from the section lift vector, which is ’uniquely defined to be 
perpendicular to the wind vector. In forcing the blades in the in-plane 
principal axis direction the effective plane of the in-plane motion may vary 
by angles of the order of 10° or more from the wind direction. The error 
caused by eliminating the section drag due to lift, therefore, will be vir- 
tually negligible in determining the in-plane aerodynamic external forces. 

In-plane motion of a single blade . - In an actual cantilever blade with 
no root motion, the flap and in-plane motions are coupled. The blade still 
possesses modes that are primarily flap motion or in-plane motion but each has 
some component of the other. 

At some feathering pitch, near zero, the two first cantilever modes, 
characterized by motions at the three-quarter radius, for example, will be 
uncoupled: pure flap and pure in-plane. The inclination of the principal 

axis will be zero. 

At zero rpm, if the hub inboard of the feathering hinge is very stiff 
so that most of the blade flexibility is outboard, then the structural prin- 
cipal axis and the flap-in-plane mode vibration principal axes follow the 
feathering pitch rotation. This is the case for the 33-ft rotor. If the hub 
inboard of the feathering hinge is flexible and the blade outboard very stiff 
then the principal axis does not rotate as the blade feathers. Most rotors 
in service approach this condition. 

At high rpm centrifugal forces come Into pjay. They produce a large 
effective stiffness resisting displacements norma] to the disk plane and a 
small effective stiffness resisting displacements in the plane of the disk, 
the latter proportional to effective in-plane pivot offset. The principal 
axis of the centrifugal stiffnesses are not affected by blade feathering 
pitch. For example, the flap in-plane principal axes for an articulated 
rotor are always parallel and perpendicular to the disk plane and independent 
of blade feathering pitch. 

Centrifugal stiffening acting in conjunction with flap and in-plane 
structural stiffnesses generates a new principal axis which rotates in pro- 
portion to feathering pitch with a factor of proportionally either greater or 
less than unity. In these analyses the factor is denoted by the symbol k . 
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Flap-in-plane inertia coupling . - The in-plane blade-hub mass equations 
of motion were written as though all mass motions took place in the disk plane 
normal to the shaft. The vertical motions equation were written as though all 
masses moved only in the vertical direction. Under these conditions no cou- 
pling between vertical motions and in-plane motions would exist. 

If a rotor possessed no precone or twist and its collective pitch and 
angle-of-attack were zero it would satisfy these conditions for small cyclic 
pitch applications. 

The 7. 5-foot rotor, studied experimentally, essentially met these require- 
ments even though it was operated at significant values of collective pitch 
and angle-of-attack. Even though the blades of the 33 foot rotor were twisted 
and preconed it is felt to have approximately met the requirements because it 
was tested only at zero angle of attack and a collective pitch, at the three- 

quarter radius, of 0 = 1.5*. Also, the blades of both rotors were stiffer 

. f 5R 

than those in use in most flying hingeless rotor helicopters. 

The vertical motions equations were forced by large, easily calculated, 
linearly independent aerodynamic and centrifugal external forces. The re- 
sulting response was stable and well behaved. 

The forcing functions for the in-plane equations, however, were nonlinear 
functions of the vertical equations external forces and response motions and 
were small in magnitude. The in-plane response was lightly damped and, near 
resonance, of large magnitude. 

High-speed gyroscope . - The high-speed gyroscope, in these analyses, was 
mounted so that It always remained parallel to the swashplate. Slop and 
elastic distortion in the linkages attaching its housing to the blade feath- 
ering horns was assumed to be zero. 

The gyroscope itself was assumed to have no vertical depth; all its mass 
was assumed to be concentrated into a flat disk. With its diametral iner- 
tia, or moment of inertia about an axis lying in the disk, denoted by 1^, 
its unforced equations of motion in stationary axis were: 
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where ft is the gyro rotation rate in radians per second 
G 

Mean aeroelastic derivatives . - If the shaft of a hingeless rotor were 
fixed so that its response motions were zero and it was then given prescribed 
steady rigid body displacements and velocities so that steady aeioiynan: 
states existed at the rotor, then the rotor would attain steady oscillating 
deformed shapes under the action of the mot ion- induced airloads ar gyroscopic 
forces. It would also develop steady mean and oscillating forces at the hub. 


The mean aeroelastic derivatives are the rates of change of each mean 
hub force component with respect to each change in shaft displacement and 
velocity. Derivatives may also be found if the change in each mean force 
component due to unit changes in each control system displacement and basic 
geometric shape such as twist and precone is given. 

Mean aeroelastic derivatives are often used in conjunction with the six 
degree of freedom rigid body equations of motion for helicopter overall sta- 
bility analyses in much the same way as wing aeroelastic derivatives are used 
in fixed wing aircraft stability analyses. 

The mean and oscillatory aeroelastic derivatives may be calculated using 
the rotor alone vertical motions equations, which in this study possess num- 
bers of degrees of freedom equal to the number of blades. The aerodynamic 
coefficients of the differential equations or the ratio of change of g r, ~er- 
alized force per unit displacement or velocity of each degree of fieedom tv st 
first be calculated. The coefficients are divided into two classes: -esp.nse 

coefficients and forcing coefficients. Response coefficients represent the 
aerodynamic forces produced by motions of the degrees of freedom. For.’. ng 
coefficients represent those aerodynamic forces produced external to +,>; dy- 
namic system by gust angle-of -attack, rotor preconing, and blade twist and by 
swashplate collective, which is not a degree of freedom in this analysis. 
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TWIST SHAPE PARAMETER 



+ twist is nose up and it is measured 

from zero at the root (g rotation). It 

should be noted that 0^ is a rate of 

change of twist angle per unit span; 

therefore 0, . = 0.R. 

tip t 



The main difference between aerodynamic derivatives of fixed wing and 
rotary wing aircraft is that the rotary wing derivatives or coefficients pos- 
sess, as well as a mean value, parts which vary periodically with time (or 
azimuth). In the case of the forcing derivatives, these merely add steady 
oscillatory components to steady-state conditions. But the periodic parts of 
the response coefficients alter the basic mathematics of the differential 
equations . 

The coefficients, analogous to the rigid wing derivatives of fixed wing 
aircraft, relate the rotor forces, moments, and generalized forces to the 
displacements and velocities in the various rotor degrees of freedom. The 
coefficients are employed in equations relative to an earth-fixed axis system, 
but are serially derived in this report from blade forces due to blade motions 
through full rotor coefficients in rotating axes to the final fixed axes values. 
Coefficients may be found for rotors with three to five blades. 
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An examination of the stationary axis coefficients showed them to consist 
of mean values, independent of rotor azimuth, and harmonic components at f re- 
frequencies of integer multiples of the number of blades times the rotor 
rotational rate. The lowest frequency multiples of the three- and four-blade 
rotors were 3P and 4p respectively, and were much larger in magnitude than 
that of the higher multiples of rotor frequency. For this reason only the 
amplitudes of the lowest frequency components, along with the mean value, 
were kept in the analyses. The phase of the harmonic components, remained 
essentially unchanged except at values of 1/V approaching zero. 

The mean aeroelastic derivatives based on rotor alone equations of mo- 
tion, including aerodynamic coefficients just described, are displayed in 
vector diagrams where changes in two moment or force components take place 
with respect to one controllable motion. For example, the variation in hub 
pitch and roll moment per unit cyclic pitch aeroelastic derivatives, with 
blade stiffness and advance ratio, is shown as a map of hub moment vector 
change . 

Oscillatory aeroelastic derivatives. - The oscillatory components of the 
aeroelastic derivatives , relative to stationary coordinates, consist primarily 
of first harmonic, bfl, frequency; where b is number of blades and Q the 
rotor rotation frequency. (Higher harmonic forces are very small and are not 
studied in this report.) 

The bO harmonic variations of hub and swashplate moment and shaft 
shear execute fixed ellipsoidal patterns in a vector diagram such as that 
used to display mean derivatives. The ellipse may be flat or circular or 
have any aspect ratio between these limits and its major axis may be oriented 
in any azimuth. For example, the shaft shear oscillatory forces produced by 
a 3-blade rotor in stationary axis may be described as the following function 
of azimuth, where Y and X are lateral and fore -aft forces respectively. 

Y = Y cos 3i|r + Y sin 3^ 
jC jS 

X = X-, cos + X 0 sin 3^ 

3c T 3s 
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This resolution of the forces may be transformed into rotating coordinates 
to yield the two constant amplitude components, 2P and 4p, of vibration which 
had combined to form the 3P ellipse in stationary axis. 

- A cos 2\|f + B sin 2i|i + C cos 4^ - D sin 4\jf 
H 

X = B cos 2\|f - A sin 2\|i + D cos 4^r + C sin 4i|f 
K 

In rotating coordinates the force vector consists of two constant ampli- 
tude parts one advancing at two per revolution (2P) and the other regressing 
at four per revolution (4P). These two force components fully describe the 
oscillatory forces, whether in rotating or stationary coordinates. 

If the two components are kept separate and transformed back into sta- 
tionary coordinates, each keeps its same advancing or regressing character but, 
of course, changes back to 3P frequency. 

Y = A cos 3^ + B sin 3t + C cos 3^ - D sin 

X = B cos 3^ - A sin 3t + D cos 3^ + C sin 3^ 

The two vector components are shown in Figure 9* 
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Figure 9- Three-Blade Rotor Shaft Shear Force Variation With Time 
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The 3P two-dimensional moment or force vibration may therefore be 
completely described in terms of its advancing and regressive components 
by specifying the amplitude and azimuthal heading, or the Y and X components, 
of each at a particular instant of time or rotor azimuth position. In this 
report the vibration components are described at rotor azimuth ^ = 0, that 
is, with the number 1 blade pointed aft along the X-axis. 

It is clear that the Y and X force components of each of the advancing 
and regressing components of the 3P vibration at i|i = 0 are given by the 
following : 


Y 


X 


2p Contribution 


Y^ - X 
= A = — — 

2p 2 


X„ + Y 
= t. = 3c 

2p B 2 


3 s 

3s 


4p Contribution 



3c 


+ X 


3s 




The complete specification of oscillatory aeroelastic derivatives of 
shaft two-dimensional forces follows the rules of shaft shear: 


a 

a 


gp’ 


lc : 





Residual forces . - The concept of residual forces facilitates the dis- 
cussion, of experimental results. In a rotor wind tunnel experiment only the 
following parameters may be controlled and varied: cyclic and collective 

nitch 0,9,0 and angle -of-attack a. It is therefore only possible to 

* c s o 

obtain experimental aeroelastic derivatives of hub forces and moments and 
blade forces with respect to these parameters. 
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In rotors with precone, blade twist, and camber, there are "residual forces" 
with the controls neutral, i.e., with 0 - 0 = 0 = y ~ 0 . Knowledge of the 

residual mean and oscillatory forces coupled with the mean and oscillatory 
aeroelastic derivatives allows the recreation of all test conditions and ex- 
trapolation and interpolation of them. In particular, cyclic pitch required 
to trim hub moments to zero may be found, and these are displayed to indicate 
the center of cyclic pitch variation during the tests. In addition, cyclic 
pitch to trim swashplate meoments or any 2P or 4p vibration component to zero 
could be determined from this data, if desired. 

It should be noted that the residual forces for the 7 • 5-ft-diameter rotor 
are zero. It should also be noted that since no systematic variation of col- 
lective pitch or angle -of-attack were made in the 33-ft rotor tests, the con- 
tributions of the small values existing during the tests are added to the 
residual forces. Because of the large cyclic pitch required to trim out the 
effects of twist and precone on the 33-ft rotor it was not possible to actually 
test at zero cyclic pitch under all rpm-forward speed conditions without over- 
loading the rotor. At these test conditions the residual forces may not be 
realistic, since they were obtained by extrapolation of a best fit plane through 
the experimental data, and at the high rotor forces at zero cyclic pitch non- 
linear effects may be important. 
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VERTICAL MOTIONS EQUATIONS 


Equations of motion are derived for a free flying rotary wing aircraft 
having three body degrees of freedom, a main rotor with an arbitrary number 
of blades (b),and a control gyro on the main rotor with two degrees of free- 
dom. The body degrees of freedom are restricted to pitch, roll, and plunge; 
the gyro degrees of freedom to pitch and roll only. Only one bending mode is 
used for each individual rotor blade. Its shape is parabolic and in many 
cases (in particular for a stoppable rotor vehicle) is a good approximation 
of the first vertical or "flapwise" blade bending mode. 

The derivation first develops the equations of motion of each single blade 
in its own axis system rotating with the rotor. The equations, in four degrees 
of freedom, contain matrices of inertia; centrifugal, structural, and aero- 
dynamic coefficients . 

Next, the independent sets of equations, one for each blade, are coupled 
together and transformed into equations of motion of the total rotor, still 
in rotating coordinates. Coefficients representing' swashplate springs and 
dampers and rotor elastic mode structural damping are then added. 

Following this, the equations in rotating coordinates are transformed into 
stationary coordinates where the gyroscope, swashplate stationary axis springs, 
and dampers and body terms are added to complete the homogenous linear ordinary 
differential equations. 

Lastly, the external forces applied to the system are calculated. These 
forcing functions, occupying the RHS of the equations, are caused by control 
motions or forces, precone, twist, and angle -of-attack. 

Single Blade Equations of Motion 

The single blade forced linear ordinary differential equations of motion 
establish the condition of dynamic equilibrium between , the external generalized 
forces and those produced by motions of the blade degrees of freedom. They 
consist of sums of the products of coefficients and motions as follows: 
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rotor precone 
rotor collective 


blade twist rate 


gust angle-of -at- 
tack factored by 
forward speed 


and the square matrices represent the following 

Bis le inertia 

Centrifugal and structural stiffness 
Aer odynami c dampi ng 

Aerodynamic stiffness 

Blade aerodynamic and centrifugal 
forcing 

Single blade inertia matrix , - The single blade inertia matrix represents 
the relationship between accelerations in the degrees of freedom and general- 
ized forces on the degrees of freedom due to the accelerations. 
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Note : 


b .m. 


B 

V 

- Tt 1 

z 

p.m. 

- - N 

e'c / 4 



_ 6 


The minus sign merely indicates 
the presence of the forces on 
the RHS of the equation in this 
expression. 


The blade generalized forces are as follows: 

b.m. Blade root bending moment at the center of rotation, ft lb. 
V Blade root shear at the center of rotation, lb. 
p.m. Pitch moment about the blade quarter chord, ft lb. 


H Blade flapping generalized force, lb. 


The 




matrix therefore is defined as follows: 


9b .m. 

9b. m. 

9b. m. 

9b .m. 

as 

a'z 

ae c/ 4 

96 

av 

av 

av 

9V 

as 

az 

ae c/ 4 

96 

9p.m. 

9p.m. 

9p.m. 

9p.m. 

as 


90 c f 4 

96 

9H 

9H 

ffl 

9H 

as 

9z 

99 c/ 4 

96 


For blades with mass centroids of sections distributed along the 
quarter chord line, mass distribution given by ^ (r) and a local pitching 

moment of inertia of the blade about the quarter chord of I , the above 


matrix becomes: 



The matrix is symmetric. 
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Single blade centrifugal and structural matrix , - The single blade 
centrifugal and structural matrix relates the generalized single blade forces 
to displacements of the single blade degrees of freedom. 





Note: Minus sign indicates terms are 

on RHS of equations. 


The matrix is therefore defined as follows: 


^b»m . 


Bb .m. 

30- /i. 


9b .m. 
96 


H 


9p.m. 

9B 


9p.m . 

9z 


9p.m. 

"cA 


9p.m. 

96 


The centrifugal part of the matrix may be formed from the inertia matrix 

2 

by factoring moments of inertia by Cl and making the second row and column 

dH 

zero. The structural part of the term — r is obtained by substituting the 

2 d ° -2 
natural frequency squared, u> 6 , for the rotational frequency squared, Cl , 

in the product with the flapping generalized mass. 
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Single blade aerodynamic matrices . - There are three single blade 
aerodynamic matrices. They relate root bending moment, shear, quarter-chord 
pitching moment, and flap generalized force, to the velocities and displace- 
ments of the degrees of freedom and to the fixed geometric shape parameters 
of the blade. 


The two response aerodynamic matrices are: 


K 1 


3F 


and 

(♦) 

3H 


and the forcing aerodynamic derivatives in matrix form are: 


9F V 


3 Fixed 
geom 


(♦) 


The two response matrices are similar in form to the centrifugal and 

For example, 


(♦) 


3F\ 


(♦) 


3T1 


; except terms 
damping matrix 

in each are functions of 
is as follows: 

azimuth. 

3b. m. ^ 
30 

3z 

3 b « m . [ . \ 
. (♦) 
30 c / 4 

3 b .m. ^ 
36 


^(♦) 

3z 

ae =/ 4 

*?<♦) 

36 

30 

3z 

^•(*) 

30°/^ 

3p.m. ^ 

*?<♦> 

30 

«<♦> 

3z 

30 c / 4 

*<♦) 

at 


The forcing matrix, on the other hand, is not square. It has four rows, 
one for each generalized force, and columns equalling the number of rotor and 
blade fixed geometry descriptive elements. The fixed shapes considered in 
this analysis have been: rotor precone B , rotor collective pitch 0 , 


blade twist rate 0 ^, and gust vertical velocity z 


g 
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In addition to aerodynamic forcing functions, the forcing matrix contains 
terms due to centrifugal force acting on rotor blade linear flapping due to 
precone and collective pitch. 


The combined aerodynamic and centrifugal forcing matrix is 


as\ 


|d Fixed 
geom 


<♦) 


o 



3b . m . 
z g 




|^-U) 

z 

g 

IT’O) 

o 

o 

3p.m. 

* 6 t 

3p.m. 

g 


3H / . \ 

*<•> 

f. <♦> 

z 


H 


as follows : 

1 sin A 0 0| 

0 0 0 0 

0 cos A 0 0 

0 0 0 0 


The elements of the aerodynamic matrices are evaluated at closely spaced 

intervals of azimuth. They represent values of root b.m., shear, quarter- 

chord pitch moment, and flap generalized force due to each of the displacements 

and velocities of the degrees of freedom. and also the fixed geometric shapes, 

and are found by integrating aerodynamic forces radially. These aerodynamic 

coefficients are functions of advance ratio ^ and tip speed dynamic pressure, 
o 

2 q. . = p(QR) as well as, of course, the rotor blade detailed geometry. 

n;ip 

Aerodynamic strip theory is used. That is, the flow relative to the 
blade is resolved into components parallel to the blade quarter-chord line 
and normal to it, section dynamic pressure is based on the normal component 


Q 

q = £ v , and the angle -of- attack of the section a is measured between 
u n 2 n n 

the normal component and the blade chord line. The effects of the radial flow 
component parallel to the blade are ignored. The section lift is given as 


df = Cj? . c* n q n cdr 

a 

where c. = section lift curve slope 

a 

c = section chord, ft 
dr = increment in radius, ft 
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Induced inflow and unsteady aerodynamics have been neglected and tip 
losses accounted for by the tip loss factor B. 

The distributions of lift due to blade motions and geometric shapes are 
integrated radially so as to yield the four generalized blade forces b.m. ? V, 
p.m. 5 and H at closely spaced intervals of azimuth. This gives the elements 
of the three matrices as functions of azimuth. 


The effects of the reverse velocity region were explicitly accounted for. 
The aerodynamic center of the blade was assumed to shift to the three-quarter 
chord point. 


In the three matrices there are only three types of section angle -of-at- 
tack; they are due to: section pitch, radial slope, and section velocity: 

section pitch 


radial slope 


section velocity 


The increments of lift in the radial direction for each type of angle-of- 
attack are as follows : 


Section pitch: 

a = 0 /), + 9 + 0, r 

n c/4 o t 


dt = c .a £ (Or + V sin il/) 2 

o n 2 v T 

ot 


This expression applies in the advancing flow region, 
velocity region the sign of the lift increment reverses. 


cdr 


In the reverse 
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R' 1 '.til slope : 


slope (radians) = 8 0 


+ 8 + 



6 


_ - (slope) V cos 
a n ~ (Qr + V sin If) 


d£ = c (-slope) V cos ^ ^ (Or + V sin t ) cdr 


a 


This expression applies in advancing flow. The sign changes in the 
r? vers. 2 velocity region. 


2 . 

Section velocity: velocity = r8+[(r) 6 + z 


• •if) 


ot. 


velocity 


n (Qr + V sin \|i ) 


d £ ~ c i (-velocity) (Qr + V sin t) cdr 

ot 

And again the sign changes in the reverse velocity region. 

The radial distribution of the lift increment dl at each azimuth posi- 
tion is then factored by the mode shapes of the four modes to yield the aero- 
dynamic derivatives. The four mode shapes are (l) linear flapping, (2) ver- 
tical displacement, (3) pitch moment arm about the quarter-chord line, and 
( 4 ) parabolic flapping. 


b .m. 


V 



dr 


dr 


p.m. 


/<=-!> SF dr 


(reverse velocity region 
only) 


H 



dr 


The above integrations at each of a large number of azimuthal positions 
become the aerodynamic derivatives for the single blade. 
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Rotor Equations of Motion 


The single blade equations would permit the calculation of the independent 
motions of each blade unrestrained at its root by the shaft or by attachment 
to the other blades through the gyroscope. In this section, the restraints, 
offered by the shaft to vertical and tilting motions and by the other blades 
(through the swashplate) to feathering motions, are applied to the three in- 
dependent blade equations causing them to be transformed into the equations 
of motion of the complete rotor in rotating coordinates. 

Following the transformation to rotor degrees of freedom in rotating 
axes, the swashplate rotating damping and rotor flapping structural damping 
terms are added. The rotor equations are then in a form to be transformed 
to stationary axes. 

Once in stationary coordinates, stationary damping, springs, and gyro- 
scope terms are added to the swashplate and the body equations of motion, 
complete with pitch, plunge, and roll inertia and aerodynamic terms, are added 
to the rotor equations of motion. 

Transformation of blade external forcing functions due to precone, twist, 
collective, and angle of attack through rotor rotating coordinates and adding 
swashplate control moments or displacements and body residual forces and mo- 
ments complete the rotor airframe equations of motion. 

Rotor equations of motion in rotating axes . - The pth single blade motions 
are represented by the vector: 



6 

P 


For a three-blade rotor, for example, there would be a corresponding 
vector for all blades 
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and it would contain 12 components, or degrees of freedom. 

It : * possible to write the equation of all three independent blade 
motions together as follows : 

fb 

311*2 


^b 

9 fixed 
geom 

^b 

9 fixed 
geom 

^b 

9 fixed 
geom 

It has been shown in the section on kinematics that blade displacements 
are related to rotor displacements in rotating coordinates by: 



Therefore i.he vector of all blade displacements is related to the rotor mo- 
tions by: 






74 



or 


and since the relationship between blade motions and rotor motions is purely 


kinematic, then: 


N 

M * 


On the other hand, the transpose of the 


matrix relates the generaliz- 


ed forces applied to the rotor degrees of freedom to the generalized forces 
on the individual blades. 
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Generalized Rotor Forces in Rotating Axes 


where M 

r 


L 


r 


T 

r 


M 


0 


r 


Rotor pitch moment ft lb 
Rotor roll moment ft lb 
Rotor thrust lb 
Swashplate pitch moment ft lb 



r 


Swashplate roll moment ft lb 


H c 


Rotor collective flapping generalized force lb 
Rotor pitch flapping generalized force lb 


H c = Rotor roll flapping generalized force lb. 
r 

The two properties of the JdJ matrix permit the three uncoupled rotor 
blade equations of motion to be transformed into equations of motion of the 
overall rotor degrees of freedom in rotating coordinates. 


76 




Rotor equations in rotating coordinates may be written more concisely 
as follows : 
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The [AF (♦)] matrix contains aerodynamic forcing functions, primarily, 
but it also includes centrifugal collective generalized forces due to precone, 
0^, and collective pitch, 0 q . 

Rotating damping, C , due to feathering friction, and rotor flapping 
structural damping, may be conveniently added at this point. They form 

a mechanical damping matrix relating rotor degree of freedom generalized 
forces to velocities of the degrees of freedom 



| generalized forces 



z 


[ 

°a] 

a 

diagonal 

matrix is defined 

as follows: 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C R 

0 
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0 

0 

0 

0 

0 

0 

C R 

0 

0 

0 

0 

0 

0 

0 

0 

2 

,J) 6 M 6 6 
o o 

0 


0 

0 

0 

0 

0 

0 


0 0 0 0 0 0 0 


2 y, 0) M 

6 6 6 Q 6 0 


0 


0000000 0 2 Y, U). M. . 

6 6 6 0 6 0 


The final equations of the rotor relative to rotating coordinates are 
as follows : 

[i] B z + [Da] f3 z + [S] 0 z - (♦)] 0 z - [A = [A.F. (♦)] f 

where 
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Transformation to stationary axes . - The transformation matrix relating 
the rotor rotating degrees of freedom, 0^, to the rotor stationary axes de- 
grees of freedom, 0^, and discussed earlier, is jT^j . 

B z = [ T z] B Z 

The transpose of the transformation matrix supplies the relationship 
between the generalized forces on the degrees of freedom in stationary coor- 
dinates to the generalized forces on the rotating axes degrees of freedom. 

jGeneralized Forces^j = j^zJ T {generalized Forces^j 

The transformation and its derivatives and transpose supply the means of 
transforming the equations of rotor motions from rotating to stationary axes. 

The derivatives of the relationships between rotating and stationary 
degrees of freedom are as follows: 

K - [*.] B z + pz] B z 

*z - [ T z] h * 2 pz] B Z + [ ¥ z] h 

The equations of motion in rotating axes may then be transformed to sta- 
tionary axes in the same fashion as equations of blade motion were transformed 
into equations of rotor motion. The full transformation process is written 
out and then the abbreviations permitted by rotational symmetry are shown: 


[ T z] T M{[ T z] h * 2 pz] B Z + pz] B Z ) 

+ PI M f PJ ^ + Pz] 8 Z 1 + Pz] 1 M [ T z] 9 Z 

- Pz] T N ( pz] 9 Z. + pz] 9 z) - Pz] T W pz] 9 Z 

= [T z 1 T [A.F.] f 

The derivatives of the transformation matrix consist of derivatives of 
the elements of the matrix. For example, JtJJ is as follows: 
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n 


sin fi t 
cos fi t 
0 

0 

0 

0 

0 

0 


- cos 0 t 0 

- sin 0 t 0 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


and similarly for 



0 

0 

0 

sin 0 t 
cos 0 t 
0 
0 
0 


0 0 0 

0 0 0 

0 0 0 

cos n t o o 

sin Q t 0 0 

0 0 0 

0 0-sin fi t 

0 0 cos fi t 


0 

0 

0 

0 

0 

0 

cos Q t 
sin Q t 


The inertia, centrifugal and damping aspects of the rotor are all rota- 
tionally symmetric and are independent of rotor azimuth position 1/ or fit. 
For this reason transformation of the matrices describing them to stationary 
axes can be simplified by employing the values of the transformation, its 
derivatives and its transpose at i|i = 0. The same result would be obtained 
if transformation matrices at any other values of azimuth had been employed. 


The transformation matrices 

at 

♦ = 

0 are 

as follows : 
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1 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 

0 

T z (o)] = [t 2 (0)] T = 

0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

1 


0 

a 

0 

0 

0 

0 

0 

cT 


1 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

r* / xi 

0 

0 

0 

0 

-1 

0 

0 

0 

It (o) = n 









l z J 
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0 

0 

1 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

-1 


0 

0 

0 

0 

0 

0 

1 

0 
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[*z <♦>] = - 


The |t z (ojj, [T z (0)] and (o)J are employed in transforming the 
inertia, centrifugal and structural, and damping matrices to stationary co- 
ordinates. The transformation and derivatives varying with azimuth |t z (i|r) j 
and |^T z (♦)] , however, must be employed in transforming the aerodynamic matrices 
to stationary coordinates since they vary with azimuth. 

Including the simplifications, the rotor equations in stationary axes 
become : 

[I] B z ♦ \z [I] [tJ + [Da]] e 2 + jjl] [Sj + [Da] |\] + [s]j 0 Z 

- [ T z ( *>] T pn][ T z <♦>] 3 Z - [ T z (*>] T [[\]f f Z <♦>] + M [V*>| 3 Z 

= [t b (t)J T [AF] f 


To the equations of motion of the rotor in stationary coordinates must 
be added terms for swashplate springs and dampers to ground. The swashplate 
rotationally symmetric spring and damper terms are as follows: 



The final rotor equations in stationary axes including the above terms 
are as follows: 

' [l] B z + [c] 0 Z + DO 0 z - [b b ] B z - [B] 0 Z = [f] f 

arid the matrices contain terms of the following kind: 
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Inertia 

Mechanical damping and gyroscopic 
Mechanical stiffness 

Aerodynamic damping, function of azimuth 

Aerodynamic stiffness, function of azimuth 
Forcing functions 

Complete Vertical Equations of Motion 

The rotor equations and forcing functions have been derived in stationary 
coordinates in the previous section. The gyroscopic terms are shown on page 
58. In this section the body inertia and aerodynamic coefficients are derived 
and combined with the rotor and gyroscopic equations to form the left-hand side 
of the equations or the homogenous differential equations. 

It should be noted that the term "body” refers to the complete non- 
rotating configuration to which the rotor is attached. It includes the body, 
nacelle, wing, and tail. 

The rotor forcing functions, the body forces and moments existing at zero 
angle-of- attack and swashplate control moments are combined to form the right- 
hand side of the complete equations or the forcing functions. 

Body Terms . - The body terms consist of inertia forces due to acceleration 
of the three body degrees of freedom pitch, roll, and plunge, and aerodynamic 
terms due to acceleration, velocity, and displacement of the degrees of free- 
dom. In addition, there are body steady forcing terms due to c.g. offset, 
aerodynamic forces at zero shaft angle-of-attack, due to gust angle-of-attack 
and jet engine thrust. 

Inertia . - The equations of motion of the rotor were derived with the 
coordinate axes assumed to pass through the center of gravity. With the c.g. 
on the shaft centerline, Z forces (lift) produced no pitch or roll moments 
about the c.g. 

It is now required to modify the equations to accept changes in c.g. posi- 
tion both fore and aft, and laterally. (The XH-51A(C) has a fairly wide lat- 
eral c.g. range due to vehicle asymmetry.) 


[1 

[c] 

[E] 

N 

[B] 

[f] 


82 



In aircraft work it is standard practice to consider the coordinate axes 
to pass through the c.g. and the aerodynamic terms are altered to suit changing 
c.g. positions. In the case of the helicopter, it appears to be simpler to 
leave the axes unchanged and modify the inertia matrix to suit the changing 
c.g. position. 

The following figure illustrates the changes: 



About the c.g. axes (X , Y*) the inertia forces are as follows: 


__ 
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Moments about the aircraft axes (X, Y) are related to moments about the 
c.g. axes as follows: 
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0 

-AX 
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By the transpose rule the c.g. axis motions must be related to the 
aircraft axes motions by 




to aircraft axes by equation (2) yields 



and this yields 

I + MAX 2 -MAX AY -MAX 

-MAXAY I , + M AY 2 MAY 

-MAX MAY M 

JXjody] ®z 

Thus, with pitch and roll moments of inertia about the c.g. axes called 
and arK * mass body (but not blades) called M and 

with the displacements of the c.g. relative to the shaft-centered aircraft 
axes called Ax and Ay? the inertia terms of the aircraft body are as shown 
above . 

Aerodynamic . - The aerodynamic terms in the airframe equations of motion 
are as follows for the unforced system: 
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= 0 





= 0 


Some small aerodynamic coupling terms may also exist but they have been 
ignored . 


Before these equations can be added to the complete rotor equations 
(those which include the forcing terms on the RHS) , the forcing and balancing 
wing -body-tail-nacelle terms must be added to the RHS of the above equations. 


It should be noted that the signs of the aerodynamic terms are such 
as to convert them from aircraft axes to helicopter axes. Therefore, the 
aerodynamic terms have been derived according to aircraft sign conventions. 



AIRCRAFT AXES 


HELICOPTER AXES 


Both are right-hand systems. 

The RHS terms consist of gravitational attractive forces, gust angle-of- 
attack forces, aerodynamic pitch, roll, and lift forces at ^ ro -^ or = the 
change in body forces due to jet engine thrust, and a force correction term 
due to rotor downwash and to correct balance error. 
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When accelerated or given a tilting velocity, the gyroscope generates 
body and swashplate moments as follows: 
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The left-hand side terms are added to the rotor body differential 
equations to account for the effects of the gyroscope. 

Complete vertical equations . - The complete equations of vertical mass 
element motion with forcing functions consist of rotor equations, body equa- 
tions, and gyroscopic equations. The complete set is as follows: 

[ T + X ho:ly + J RR body + B Z + [ C + J R body + C g] B Z 

+ [ E + J body] B Z - [ B R ( *>] B Z - [ B <*>] B Z 

- [f (*)] (f| + (F body ) + |cm| 

Swashplate control moments jc.m.| have been added. They are applied in a 
vector with all elements zero but the Mg and M ^ swashplate moments. Ex- 
ternal control moments may be applied to the free swashplate through them for 
,f closed loop" operation. 

The equations represent the free flight of a feathering feedback gyro- 
scope- stabilized hingeless rotor compound helicopter. 
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IN-FLANE MOTIONS EQUATIONS 


Equations of motion of the rotor blade - hub mass system in the plane of 
the disk or the plane normal to the shaft are derived. The system is free in 
the disk plane, or constrained by springs to ground. The hub mass may translate 
longitudinally and laterally and it may rotate relative to the rotating axes. 

Each blade is assumed attached to the hub at a pivot point located at a fraction 
of the tip radius from the shaft. The distance depends on the blade elastic 
dynamic in-plane mode shape. The blade in-plane motion about the pivot is as- 
sumed to be rigid. 

The number of degrees of freedom in the system is three for the body freedoms 
and one for each blade; i.e., 6 for the three-blade rotor and 7 for one with 
four blades. The rotor elastic degrees of freedom are collective (or all blades 
together to give a shaft torque) lateral, longitudinal, and the reactionless or 
differential collective or scissors mode. 

The derivation first develops the equations of motion of a single blade in 
its axis system rotating with the rotor. The equations, in four degrees of 
freedom, contain matrices of the following coefficients: inertia or acceleration 

terms, coriolis or velocity terms, centrifugal and structural or displacement 
terms, and aerodynamic velocity and displacement terms. The response aerodynamic 
terms are assumed to depend on the blade section drag coefficient at zero lift 
and are very weak. The periodic parts may be ignored without significantly af- 
fecting the blade response. This assumption reduces the equations of blade 
motion to linear ordinary differential equations with constant coefficients. 

The individual blade degrees of freedom are then constrained to take up 
only those motions permitted by the overall rotor degrees of freedom. The re- 
lationship has been discussed in the section on kinematics. The transpose of 
the kinematic relationship matrix is then allowed to act on the individual blade 
generalized forces so as to gather them into generalized forces on the overall 
rotor degrees of freedom. In this way the three or four sets of individual blade 
equations of motion are transformed into a single set of equations of motion of 
the whole rotor in rotating axes. Terms are then added to account for the mass 
and moment of inertia of the hub and blade structural damping. 
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Following this, the equations are transformed into stationary coordinates 
where non-rotating body mass and springs to ground (if tunnel mounted) are added 
to complete the homogeneous or left-hand side of the differential equations. 

The in-plane equations forcing functions or forces external to the homogeneous 
set are next derived. The major portion of these forces is produced by forces on, 
and motions of the degrees of freedom of the vertical motion equations. The re- 
maining small portion comes from the azimuthal variation of blade section drag. 

The major forcing functions are divided into two parts: vertical motions 

induced and vertical aerodynamics induced types. 

The in-plane forces caused by blade vertical motions come about because in 
actual fact blade element masses also move slightly in the radial direction as 
they take up vertical deflections. These small radial displacements and veloc- 
ities cause in-plane centrifugal and coriolis forces - which are considered to 
be external to the in-plane equations. 

Essentially vertical section air forces actually have small components in 
the plane of the disk and in the direction of the principal axis of blade lead- 
lag motions. These components are carefully calculated to produce the other 
major portion of the in-plane forcing function. 

Single Blade Equations of Motion 

The single blade freedoms are discussed in the section on in-plane motions 
of mass elements. The vector of dispacement of the degrees of freedom of the 
pth blade is: 



The homogenous differential equations of the in-plane motions of a single 
blade are as follows: 
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I 1 + 

COR 

f| + 

CF 

p l 

P 


P 







= 0 


where the square matrices represent the following: 


COR 

P 


CF 



A 

P 


Blade inertia 
Coriolis coefficients 
Centrifugal and structural stiffness 
Blade aerodynamic damping 

Blade aerodynamic stiffness 


Forcing functions are treated in a later section. 


Blade inertia matrix . - The blade inertia matrix represents the relation- 
ship between accelerations in the degrees of freedom and generalized forces 
on the degrees of freedom due to the acceleration. 


f n -s 



I 

p 


U, 



Note: Minus shows terms 

on R.H.S. 


The blade generalized forces are: 

n.f . Normal force acting on blade, or pivot shear, lb 

a.f. Axial force acting on blade, or pivot tension, lb 

s.t. Shaft torque, or root in-plane bending moment, ft lb 

p.t. Pivot torque, or in-plane moment at pivot, ft lb 
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The I matrix represents the rate of change of generalized force per 


unit acceleration. For example the (l,l) element is 
matrix is as follows: 


dn.f . 


. The complete 


M, r 
b c.g. 


% < r =. E .- e ) 


lb r c.g. 


shaft 


T - e M r 

D D C .g. 


v r - e, 

c.g. 


shaft 


eM. r x b . , 

D C.g. piVOt 


The Ip matrix is symmetrical. 


Blade mass = 


J dr 


Note: e = pivot offset 

in feet. 


Center of gravity of blade measured 
from center of rotation 


shaft 


= Moment of inertia of the blade 
about shaft centerline 


/ R 

dm 

dr 


I = Moment of inertia of the blade 

pivot about the pivot 

Blade coriolis matrix . - The coriolis matrix relates blade generalized 
forces to velocities of the blade degrees of freedom. 


Note: Negative indicates 

R.H.S. 


COR = 2Q 
P 


-M, r 
b c.g. 

-“b (r c.R.' e) 


M, r 

D C.i 


The coriolis matrix is anti-symmetric. 


,r - e 

c.g. 
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Blade centrifugal and structural matrix , - The centrifugal matrix relates 
generalized forces to displacements of the degrees of freedom and includes 
blade structural stiffness effects as well as rotating mass effects. 


n.f . 


a.f . 


s . t . 


p.t. 


CF 


CF 


d n.f. 
du 

d a.f. 
du 

d s . t . 
du 

d p.t. 
du 

0 


A, 


CF 


c.g. 


f 

i 


u i 

d n.f. 
dv 

d a.f. 
dv 

d s . t . 
dv 

d p.t. 
dv 


n \ 


n \ < r c. K .- e) 


pivot 


d n.f. 
du 

d a.f. 
du 

d s .t . 
du 

d p.t. 
du 

2 

fi^r 

0 
0 
0 

+ er 


c.g. 


Note: Negative 

indicates 

R.H.S. 


d n.f. 

dC 

d a.f. 
d£ 

d s.t. 

dC 

d p.t. 

dC 


Matrix for 
one blade 
only. 


0 M (r - e) 
d c.g. 


0 

0 


n 


c.g. 


uu. 


shaft 
J b . 


= 0 pivot 


The 


terms are all self-evident except perhaps ‘ 


and 


d n.f. 
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Their derivations are as follows: 



There are two types of forces produced by Q. The force normal to the 
undeformed position of the blade 

n 

d^ n.f. n 2 dm 

-r? =v = r w — sm 1 

d^dr dr 

and the force normal to the deformed blade needed to obtain pivot torque 

d 2 ( n . f . ) J n 2 dm . ( , • \ 

d£dr dr J 

Now it may be noted that £(r-e)=ir so that i ~ Q (l - and 
£ * i = — £ . Employing these two relationships and for Q and i small 


2 

d n.f. 
dCdr 


r Q 2 p (1 - *) 

dr r 


n 2 dm / x 

= 0 v— (r - e) 

dr v ' 


and d ( n . f . ) A n 2 dm e 

— L d = r n — - 

d£dr dr r 

2 dm 

= - e Cl — and the minus sign denotes a retarding force to 
forward rotation about the pivot. 

The two derivatives with respect to Q are therefore 

f R 

d n.f. _ n2 I / v dmdr n 2 , x 

__ -If |(r-e) - =0 ^ ( r c .g. - e) 
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d p * t • 

dC ~ ' w 


r 

a. 


or 


= n 


r 


/ \ dm 

e (r - e) — dr 
dr 


dm 


+ Q 


= - n 


(r - e) (r-e) ^ dr 


Q 


/ r ( r " e ) dr 


- i. 


+ e 


pivot 


shaft 


shaft 


N* 


pivot 


*b 
+ M b 


c.g. 


e r 


c.g. 


dH 

It should be noted that the above is not exactly analogous to the 
term of the vertical blade motion equations. It is, therefore, not possible 
to include structural stiffness in exactly the same way. In the in-plane 
formulation the pivot spring is always the same: ft lb/radian. Its value 

therefore, will not change with rpm. At zero rpm or Q = 0, the in-plane fre 
quency is given by 


so that 



d p.t. 
d£ 


shaft 


pivot 




e r 


c.g. 



Q = Oj 


pivot 


Blade aerodynamic matrices . - For the purposes of this study all aerody- 
namic forces caused by in-plane displacements of the blade will be assumed 
to be zero. 


A 


= 0 


The only in-plane aerodynamic force which is not related to the essen- 
tially vertical aerodynamic forces is the c^ or section drag at zero lift. 

o 

Abbot and Von Doenhoff indicate the c, of the IJACA 0012 airfoil to be .006. 

d 

o 

For purposes of this study, to account for roughness and to partially compen- 


sate the missing 


x cd 


_ an effective value of c, 

dc 2 o 

figure suggested 1 by N.B. Gorenberg. ) 


= .011 will be used. (This 
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It is now possible to approximate the aerodynamic forces due to blade 

in-plane motions. It will be assumed that the only forces on the blades due 

to blade motions will be due to c, . 

d 

o 

The blade aerodynamic damping motion 

page. 

Calculation of the elements of the blade aerodynamic damping matrix is 
straight forward. The derivation of the - P ‘ — * term will be shown as an 

9 ; 

example . 

Aerodynamic drag (or negative force) per unit radius due to unit Q 
velocity is as follows: 


is as shown on the following 



Due to minimum drag coefficient c^ the rate of change of normal force 
per unit radius is: 

d = - c^ — p (Or + V sin if + £ (r-e ) ) 2 c + c^ ^ P (Or + V sin i|f ) 2 c 

o o 

= - c 73 - p (Or + V sin \|r)^ + 2(Qr + V sin t) f(r-e) 


* 2 2 1 2 
+ f (r-e) c + c^ £ P (0 r + V sin ^) c 

o 


*2 2 

but £ (r-e) is negligibly small compared to the other terms. 


d n.f . 
dr 


- c^ ^ P 2 £(0r + V sin i|r) (r-e) 


d n.f. = - c pc (r-e) (Or + V sin t) 
d £ dr o 
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Single blade aerodynamic damping matrix: 
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(r-e) (r-e) (Or + V sin i| r) dr 


and the derivative 



• t . 


becomes : 


d p.t. 







V sin ij / 


In-plane Equations of the Rotor 

Once the single blade equations of motion are available, they may be 
assembled into a single uncoupled matrix. First, form the vector of all 
single blade displacements: 



for a three-bladed rotor 


The equations of motion of the three uncoupled blades become: 



The displacement derivatives assumed to be zero were not included. 


The equation may now be transformed into rotor coordinates employing the 
relationship between blade and rotor coordinates derived in the kinematics 
section, namely: 

98 



71 = D 0 

xy p xy p xy 

or expand to three blades 



and since the relationship is kinematic 

f| = Id 0 

xy | xy xy 


and 


Tj = Id I e 

xy | xy| xy 

and employing the transpose of the D matrix for the purpose of 

[ xyj 

gathering up single blade generalized forces and converting them into gener- 
alized forces applied to the rotor degrees of freedom. The transpose relation- 
ship is as follows: 
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where the generalized rotor forces in rotating coordinates are as follows: 



Lateral force, lb 
Longitudinal force, lb 
Yawing force, ft lb 

Collective lead-lag generalized force, lb 
Lateral lead-lag generalized force, lb 
Longitudinal lead-lag generalized force, lb 


These two properties of the D^J transformation matrix allow the three 
uncoupled sets of single blade equations to be transformed into rotor equations 
of motion as follows : 



With the equations written in this form it then becomes possible to in- 
clude rotating hub mass terms and blade structural damping terms. The terms 
for the rotating hub motions and rotor elastic mode structural damping are as 
follows and may be combined directly with the rotor equations in rotating 
coordinates : 
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The hub inertia terms and the rotor structural damping terms are defined 
as follows : 


“r 

J R 

V 


UU 


NL r , M , M r r 
£o£o ? QyQy QxQx 


Hub mass, slugs 

2 

Hub polar moment of inertia, slugs ft 
Fraction of critical damping 

Natural frequency of mode at rpm, radians/second 
Generalized masses of rotor collective, lateral and 
longitudinal lead-lag degrees of freedom. 


And the combined equations may be written more concisely as: 

[I I 3 + [COR I 6 + [CF 0 [A_ 1 0 =0 

| xy] xy xy| xy xy xy “ | 1* | xy 

In a similar way to the transformation from blade to rotor coordinates 
the equations may be transformed to stationary axes by employing the trans- 
formation relation between rotating and stationary axes and remembering that 
it is a function of time so that; 
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xyl xy 


3 + t 0 


xy] XY I xy] XY 


xy XY 


9 + 2 T 9 + T 9 


The transpose which relates the generalized forces in stationary axes to 
those in rotating is as follows: 

GF = It | t GF 
XY xy xy 


The 

rotor 

and 

hub 

equations 

in 

stationary axes 

become : 

It | t 

li | 

T 

8 


T 

T 

[COR 

It I + 

2 It 1 t 

l w) 

|xy| 

*y] 


XY 

xy] 


xy 

1 xy| 

1 *yj 


1 I T 8 

xy xy xy XY 


+ T CF T + T 
xy xy xy : 


T I B 


COR T + T 


xyl | xy| 


xyl xy xy "XY 


xy XY xy 


T 0=0 
xy XY 


And with the appropriate simplifications due to rotational symmetry the 
equations become, through the use of It I, [t I and I T 1 at f = 0, 


xy | xy| 


xy xy 


10+ COR +21 


xy| xy XY 


+ CF + COR T + I 


xy] xy] | xy] I xy| XY 


xy x y 


A- T 0 =0 

1 xy| XY 
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Body terms . - Terms for the body mass and spring restraint to the tunnel 
are combined with the rotor equations. 

The shaft may translate laterally f, y" and longitudinally "x" or rotate 
"u." No other degrees of freedom involve motions of the shaft. Masses of 

shaft or body attached to it therefore can cause inertia forces only in acting 

on those degrees of freedom, the forces Y, X, N. 

In the case of the N force the shaft is presumed to be torsionally unre- 
strained by the transmission and engine - as though it were fluid coupled. 

The engine applies a mean torque, only, to the shaft and free torsional oscil- 
lations of the shaft are allowed. 

The Y and X forces due to accelerations x and y will be due to the 

mass and moments of inertia of the helicopter body. 

It is assumed that effective non-rotating masses exist in the plane of 
the disk. One resists lateral accelerations and the other longitudinal. Their 
magnitudes are different because they are based on rolling and pitching mo- 
ments of inertia of the body as well as the body mass. 

Springs restraining the shaft against lateral and longitudinal motions 
are also employed. The terms for body or shaft inertia and springs may be 
combined with the rotor terms and are as follows: 



With the rotor (and rotating hub) removed these equations govern the roll 
and pitch or, more accurately, the lateral and longitudinal natural frequencies 
of the model in the wind tunnel. 


The homogeneous equation . - Including the body mass and spring terms with 
the rotor terms yields the full set of homogeneous in-plane equations. 


[ XY S XY 




ft 


XY 


! C0 R h| ®XY + bci 
- M SjCY = ° 
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In-Plane Forcing Functions 

The in-plane forcing functions are produced by vertical motions of the 
blades and through components of the vertical response and forcing aerodynamics 
vector. Before discussing the aerodynamic forces causing blade in-plane motions 
it should be noted that forcing aerodynamics in the in-plane direction are ap- 
plied to the blades due to collective pitch, angle-of-attack, cyclic pitch, 
precone, twist, and the flapping deflections and velocities accompanying them. 
These in-plane components of the essentially vertical aerodynamic forces 
depend on two considerations: 

1. Magnitude of the essentially vertical force 

2. Angle between the lift resultant vector and the effective plane of 
the blade in-plane response. 

Because the in-plane velocities are relatively small, of the order of 
10 ft/sec or so, they will not significantly alter the magnitude or direction 
of the forcing aerodynamics. 

The vertical motions equations forcing and response aerodynamics combine to 
yield net local angles-of-attack relative to the flow component normal to 
the blade leading edge. The local lift force at the blade element, if as- 
sumed to be normal to the relative wind at that section (ignoring the drag 
contribution to the aerodynamic resultant force) has a component in the di- 
rection of the in-plane response of that blade -- calculated in these analyses 
in the vicinity of the blade section at the 3 /^ radius. 

These in-plane (of the mode) aerodynamic forces along each blade may be 
factored by the rotor mode shape for that blade, integrated and combined with 
the similar generalized force contributions from the other blades to yield 
the in-plane rotor mode generalized force. The generalized force varies with 
time as the rotor turns through the azimuth and contributes to the forcing of 
that in-plane mode. The same forces can also be factored by each of the other 
rotor mode shapes to obtain the generalized forces in all rotor in-plane modes. 
Body or shaft modes, or degrees of freedom, however, displace in a direction 
normal to the shaft or vertical axis so their generalized forces must be cal- 
culated relative to the disk plane rather than relative to the principal axis 
of the blade -alone in-plane motion. 



Before leaving the aerodynamic in-plane forcing functions it should be 

noted that there is an aerodynamic contribution that is independent of "z M 

motion and forces. It is due to the drag coefficient of the section at zero 

lift c^ . Its generalized force is considered separately, 
o 

In addition to the aerodynamic contribution to the in-plane forcing func- 
tions due to rotor response, there are inertia or mass-induced forces due to 
blade vertical flapping motions. These depend only on the built-in precone 
angle and the sum of the contributions to the flapping motion of each 

blade from the rotor flapping degrees of freedom. The in-plane forces are 
due solely to displacements, velocities, and accelerations in the radial di- 
rection of the elements of blade mass distribution. The inertia forces can 
be factored by in-plane mode shapes to obtain generalized in-plane forces 
which are then combined with the aerodynamic external generalized forces to 
yield the net forcing function to the in-plane degrees of freedom. 

In-plane forcing due to vertical motions . - The vertical motions degrees 
of freedom are assumed to provide displacements parallel to the shaft only, 
in framing the vertical equations of motion. In actual fact, however, elements 
of blade mass take up small motions in the radial direction when the blade 
bends vertically in its parabolic mode shape, especially when in the presence 
of built-in blade precone angle. Figure 10 shows the small inboard (negative) 
radial displacement f, v ,r that accompanies the vertical displacement f *6 M in 
conjunction with built-in precone 



position of blade above plane = g Q r + 6^ (r) 

Figure 10. Blade Radial Displacement Due to Flap Bending 
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The displacement, velocity, and acceleration of a lump of mass at station 
"r" due to flapping displacement are as follows: 



The above is the motion of station lf r M due to flapping deflection, veloc- 
ity, and acceleration, in rotating blade axes, and is non-linear. 

The non-linearity causes no great difficulty in this analysis since the 
mass element motions are used only for producing external forces to the linear 
in-plane equations of motion. 


With the radial displacement, velocity, and acceleration of elements of 
mass of the blades known in terms of tip displacement, velocity, and acceler- 
ation it is necessary to find the motions of all blade tips in terms of mo- 
tions of the rotor degrees of freedom in stationary axes. Then with blade 
radial motions known in terms of motion of stationary axes degrees of freedom, 


it is only necessary to find the single blade generalized in-plane forces due 

to single blade radial motions and collect them by means of the |d^_| T matrix 

f IT 

and transform these to stationary axes by the 


xyj 


matrix. The forces may 
then be used in conjunction with the aerodynamic forcing functions to find 
the in-plane response. 


The motions of the blade tips are determined from the vertical motions 
response degrees of freedom by employing the transformation from stationary 
to rotating axes and from rotor coordinates to blade coordinates, as follows: 


T1 


z 




and since (3 is known as a function of time, then f| and T| may be ob- 
Z z z 

tained by differentiations of the above expression. From T| , f| s and T| the 

z z z 

required blade tip displacements, velocities, and accelerations may be extracted. 
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h 

The p single blade generalized forces due to radial motions may be 


written in matrix form as: 


8 n j o 

n T I. b 


n 2 e l 2 flpj. 

3 — J 


R 

b 0 


3 ? 


P I 4 4 

P o D — R 


3 R J b 


where I, 


r R 

I 2 dm 

~J r dr 

e 

f R 3 dm 

We - 

f 4dm 

l r ^ 


4 dm , 
r — dr 
dr 


2t Vb 8 nK 6 

r 2 3jr 

2QB o flK 

J, | 

R b 3-TT 


'b e 


2 , x dm 
r (r-e) — dr 
dr 


3 / \ dm 

r (r-ej — dr 
x dr 


6 6+6 
P P P 


and 6 is the tip deflection of the p blade 5 it is a function of azimuth 
\|i cr time Ob. 

The generalized forces on the rotor degrees of freedom in stationary axes 
are then found by f n.f._ } 


GF 

= 

m T 
T 

Id 

XY 


xy 

| xy 



p . t • ^ 
n.f 


for a 3-bladed rotor. 
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In-plane forcing due to vertical aerodynamic forces . - Three rotor forces 
contribute to the in-plane aerodynamic forcing functions, lift or thrust, 
pitch, and roll hub moments. The blade forces contributing to these three 
rotor forces are assumed to be concentrated at the blade 3/4 radius. The lift 
is assumed to be uniform around the azimuth, the pitch and roll moments smooth 
sinusoidal distributions as shown in Figure 11. 

The assumption that rotor pitch and roll moments are applied sinusoidally 
is a fairly good assumption. Should the forces in fact be applied at a dif- 
ferent radial station, the effect on the in-plane forces would be small. For 
example, if applied at blade tips, the blade forces would be smaller due to 
increased radius but the in-plane moment would not change significantly be- 
cause the in-plane moment arm increases to compensate. Should the forces be 
applied other than sinusoidally, it is unlikely that there would be a signifi- 
cant change in the IP or 2P components of in-plane force. 

The assumption that the lift is applied at the 3/4 radius is fairly good 
at low advance ratio. At higher advance ratio -it could be a poor assumption. 
If the lift were applied, say, at the blade tips, the in-plane forces would 
be much larger than if the lift were applied at, say, the $ 0 % station. Fur- 
thermore, the lift could be applied in a 2P distribution, or humps at opposite 
sides of the disk; or it could wander radially as it traversed the azimuth 
and still be a steady lift. 

These assumptions are used here because they simplify the analysis and 
certain evidence suggests that, in fact, they are not bad for many cases. 
Aerodynamic forcing is generally a smaller contribution to in-plane general- 
ized forces than are the coriolis forces and, therefore, can afford to be 
less precise. 


These assumptions allow the loading per blade to be written immediately 
in terms of the Lift T, Pitch Moment M, and Roll Moment L. The lift on the 
pth blade becomes: 


L 

P 


T 

b 


2M cos (\|f + \|r ) + 2L 

.75Rb P -75Rb 


sin 


(♦ + 1 ) 


where p is the pth blade numbered clockwise viewed from above. 
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With the blade loads specified at each azimuth point, it is only necessary 
to resolve them into components normal to and parallel to the inclined in-plane 
mode principal axis to provide the generalized forces for the overall rotor 
lead-lag modes and to the plane normal to the shaft for overall rotor rigid 
body modes. For the former, the in-plane axis will be assumed to be inclined 
proportional to the feathering displacement with the constant of proportionality 
dependent on the rotor rpm. 



The lift at a section is produced by two classes of angle of attack: 

Class I - due to rotation of the section relative to the disk plane or shaft, 
and Class II - due to precone, flap deflection, flap velocity, and pitch, plunge 
and roll rate of the disk. Only Class II angle s-of-attack produce an inclina- 
tion of the lift vector relative to the shaft axis. 


Under the assumption that the c^ is accounted for separately and the 
section drag due to section lift is negligible, the lift (and resultant) 
vector is normal to the flow at infinity relative to the sections. 

For Class I angles of attack the lift, therefore, is directed parallel 
to the shaft axes. Class I angles of attack are: 


0 , e_ , 

o lc’ 


0_ , and 9 n 
Is 7 t 


With the application of Class II angle s-of-attack, even in the presence 
of Class I lift, the total lift vector is inclined relative to the shaft axis. 
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The inclination is equal to the sum of the Class II angle contributions at 
the section. These angles are due to: 


^0, 6^, $ in stationary axes. 


Contributions to Class II angle s— of— attack from rotor flapping velocities 
in stationary axes 6 q , 6^, 6^ are neglected in these analyses because only 
steady cases are considered. Aircraft pitch attitude effects are included 
with the z/V parameter. 


The Class II angle-of-attack at station ,f r ,f due to the above is as fol- 
lows : 



The angle-of-attack is the angle between the wind remote from the section (at 
infinity) and the chord line regardless of the direction of the wind, advancing 
or reverse flow. This results in one equation for each contribution which ap- 
plies in both advancing and reverse velocity regions. The Class II part is 
that part between the wind at infinity and the plane normal to shaft "cr . M 

The basic blade motions and fixed geometry shape that contribute to the 
Class II angle of attack, obtained from the stationary axes degree of freedom 
motions, are as follows and depend on the azimuth location of the blade as 
follows : 
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6 : 

3 : 


*n = 


V cos jf 
Or 4- V sin \|r 


°^n ® Or + V sin i|r 


z : 


*n = 


Or + V sin 


6 : 



The net value of a due to all Class II contributors then is obtained 

n 

by factoring the blade motion column matrix by the Class II angle of attack 
row matrix. 

Class II angle-of-attack at station "r ,f on the pth blade: 

I 


a 


= - l 

n Or + V sin t ) 

p V 


I * 1 

i 1 ■ 


V cos (t + ^ )i V cos ( t + ♦ ) ir* 1* \ 

P i P l | *\ R 

n i ■ i 


In determining the a the value of radius "r" to use in these 

np 

analyses is r = . 75R (by basic assumption). 


It is now possible to calculate the in-plane aerodynamic forcing func- 
tion. On the pth blade with aerodynamic "normal" forces indicated by W L ff 
and the angle -of- attack relative to the vertical axis, or shaft axis called 
a and the inclination of the in-plane principal axis given by *0^, the 
in-plane force at the 3A blade radius is given by: 


In-plane force 


F. = L sin 

ip p 

P 


in fa + « 


is shown resolved into the in-plane principal axes in Figure 12. 
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Figure 12. Resolution of Section Lift into Blade In-Plane Principal Axis 

It is assumed that a + k 9_ remains less than 10°, so that 

n f 

P 

sin (« n + «9 f ) » (« n * «0 f ), and F = L (a + «9 ). 

P P P P ^ 

The last piece of information needed to complete the calculation of the 

is the value of 8^, . It is obtained from the Tl z matrix by abstracting 

1P P 

the 6 , the feathering pitch, and adding it to the collective pitch and 

blade twist at the three-quarter radius 9^. 75R‘ so that: 

0 =0 + 0 + 0,(.75R) so that: 

f pot 

P 

In summary, the expression for F. is as follows: 

ip p 


F iP = L P K + K0 f > 

p p 
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where 


l p = 


- 1 


2M 


• 75Rb 


COS 


(♦ + 


V + 


2L 


-75Rb 


sin (i ^ 


V 


0L 


”p 


Or + V sin ( + tjF ) 


i2r 


l 

V cos ('jf+ili ) i —7) V cos ('|r + 'l' )|r 
P i r 2 P | 


i'tefl 


k - In-plane principal axis factor 

e = e + e + ,75 r e 

f p o t 

p 


6 

P 


With the in-plane forces on each blade at the ’i/b radius evaluated at 
closely spaced azimuth positions, it is then necessary to organize these into 
generalized forces on the blades and collect these into generalized forces on 
the in-plane degrees of freedom in stationary axes. 


It should be noted, however, that in the foregoing discussion the gener- 
alized forces on the blade in-plane degrees of freedom were found. This free- 
dom was inclined to the plane normal to the shaft by the angle k and the 

P 

component of aerodynamic force in that direction was employed in calculating 
the blade generalized force. 


For the rigid body degrees of freedom y, x and v, however, the blade 

force component contribution to the generalized forces must be in the plane of 

the disk normal to the shaft. It is therefore necessary to calculate two F. 

IP 

one for k = 0 and one for k. 


This can be done as follows: 
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and for the whole rotor these become: 



In-plane forcing due to minimum drag coefficient . - The in-plane forcing 

function due to blade section minimums drag coefficient c^ depends on the 

o 

normal component of dynamic pressure q^ and the blade chord c. For the pth 
blade, 



c , c q 
d n 


c d ° 2 
o 


^Qr + V sin (\|f + ^ )j 


2 


n.f . 



d n.f. 


dr 


dr 


a.f . = 0 
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s .t . 


d n.f, r dr 
dr 



Form the IGF,. , | as a column 

( blade \ 



The forcing functions on the stationary axes degrees of freedom are available 
from 
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Complete In-Plane Equations of Motion 


The complete set of in-plane differential equations include the 
homogeneous equations in stationary axes and the forcing functions produced 
by vertical motions, vertical equations aerodynamic forces, and section drag 
coefficient . 


Sty S xy + l COR wl »v, + cf 




XY 


XY 


8 . 


XY 


XY( vertical motions 


a r jcy b xy ' Itcyj b xy j GF 

+ IgF I + IgF 1 

I XYj vertical aero | XYldrag coefficient. 
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SOLUTION OF VERTICAL MOTIONS EQUATIONS 


Linear ordinary differential equations of the vertical motions of the 
rotor-gyroscope-airframe system have been derived in the foregoing section 
for three-and four-blade rotors* Those coefficients of the equations which 
depend on aerodynamics vary harmonically with azimuth at frequencies based 
in the rotation rate of the rotor. The external forcing functions also vary 
harmonically with azimuth. 

In this section the equations are solved for the detailed harmonic or 
vibratory response of the degrees of freedom. From these and the external 
forces producing them, the shaft- and swashplate-transmitted vibration forces, 
the azimuthal variations of tip path displacement and, in the case of the 
free swashplate, swashplate harmonic tilt displacement wobble are determined. 
Calculated and experimental results are compared where possible. 

General Discussion 

The solution of systems of linear ordinary differential equations of 
motion possessing coefficients which vary harmonically with time is somewhat 
more complicated than solving equations with constant coefficients. It is 
the purpose of this section of the report to discuss the logic and 
methods which permit the solution of such equations for their steady 
oscillatory response to harmonically varying forcing functions. Such solu- 
tions yield shaft -transmitted vibratory loads and blade load histories. 

The basic logic will be illustrated by treating the simplest represent- 
ative system of the type: an undamped single degree of freedom with period- 

ically varying spring stiffness and forced at the period of the spring 
stiffness oscillation. 

Consider the equation: 

w 4- (A 4* B cos t) W = cos t 
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where A and B are constants and w is the dependent variable. The solution 
w must contain oscillations associated with the period of the oscillator and 
the period of the coefficients. Since both are the same, w may be represent- 
ed by a Fourier series, if subharmonics are excluded. (later it will be 
shown that subharmonics cannot exist in a true steady state for these 
equations ) . 

Let 

w = a + a_ sin t + a_ cos t + a^ sin 2t + a„ cos 2t + .... 

o Is lc 2s 2c 

then 

w = a cos t - a_ sin t + 2a^ cos 2t - 2a^ sin 2t 

Is lc 2s 2c 

and 


w = - a sin t - a_ cos t - 4a_ sin 2t - 4a^ cos 2t ... 

Is lc 2s 2c 

substituting w and its derivatives in the equation yields: 

-a sin t - a cos t - 4a 0 sin 2t - 4a_ cos 2t 

Is lc 2s 2c 

+Aa 0 + Aa sin t + Aa cos t + Aa sin 2t + Aa cos 2t .... 
is lc 2s 2c 

2 

+Ba Q cos t + Ba cos t sin t + Ba cos t + Ba cos t sin 2t 
is lc 2s 


+Ba 2c 

cos t 

cos 2t 

+ Ba_ 

’ 3s 

COS t 

sin 3t + Ba cos t cos 3t 

3c 


+Ba 1.3 

cos t 

sin 4t 

+ Ba 4c 

cos t 

COS 4t 1 

COS t 


It is important to note that each term, in which products of trigono- 
metric functions are found, may be simplified to the sum at two terms of 
single functions by trigonometric identities so that the above equations may 
be simplified to: 

- sin t a - cos t a, - 4 sin 2t a. - 4 cos 2t a^ 

Is lc 2s 2c 

Aa + A sin t a + A cos t a. + A sin 2t a„ + A cos 2t a„ 

o Is lc 2s 2c 


120 



B cos t a Q + B(^ sin 2t) 
+ ^ sin t) a^ + B(J cos 
+ B(^- cos 4t + |r cos 2t) 
+ B(^ cos 5t + ^ cos 3t ) 


a + B ("2 + \ cos 2t) a lc + B(^ sin 3t 

3t + ^ cos t) a 2c + B(| sin 4t + ^ sin 2t) a 


a^^ + B(|- sin 5t + ^ sin 3t) a^ s 


a t — cos t 

4c 


All terms in common frequencies are arranged in separate equations in 
order to permit a solution for the coefficients. Equating coefficients of 
common frequency: 

(a o } (a is } (a ic } (a 2s } (a 2c } (a 3s } (a 3c } 


Aa 


k B a. 


lc 


(A-l)sin t a \ B sin t a 0 

IS rS 


B cos t a (A-l)cos t a. i B cos t 

o lc Pc 


B sin Pt a, 


Is 


(A-4) sin Pt 


b B sin Pt a, 


i B cos Pt a_ 

? lc 


3s 


(A-4) cos Pt a. ir B cos Pt 

v ' Pc ^ 3c 


= 0 
= 0 
= cos t 
- 0 
= 0 


In matrix form this becomes: 


A 

0 

B/2 

0 

0 

(A-l) 

0 

B/2 

B 

0 

(A-l) 

0 

0 

B/2 

0 

(A-4) 

0 

0 

B/2 

0 


0 

0 


a 

o 


’ o' 

0 

0 


a n 

Is 


0 

B/2 

0 

< 

a ic 

► _ < 

1 

0 

B/2 


a 2s 


0 

(A-4) 

0 


a 2c 

1 

0 


It is obvious that the matrix could be extended to any desired number 
of harmonics by analytic continuation, as in the following example, to four 
harmonics. 
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~A/2 

0 

B/2 

0 

0 

0 

0 

0 ~~ 


^ 2a 

o 


0 

0 

(A-l) 

0 

B/2 

0 

0 

0 

0 


a T 

Is 


0 

B/2 

0 

(A-l) 

0 

B/2 

0 

0 

0 


a lc 


1 

0 

B/2 

0 

(A-4) 

0 

B/2 

0 

0 


a 2s 


0 

V 

0 

0 

B/2 

0 

(A-4) 

0 

B/2 

0 


a 

2c 


0 

0 

0 

0 

B/2 

0 

(A-9) 

0 

B/2 


a 0 

3s 


0 

0 

0 

0 

0 

B/2 

0 

(A-9) 

0 


a 3c 


0 

0 

0 

0 

0 

0 

B/2 

0 

(A- 16 ) 


A*. 


0 

4 J 


It should be noted that the matrix is symmetric and also that the sine 
and cosine terms are not coupled to one another, thereby allowing the equa- 
tions to be rewritten as follows: 


A/2 B/2 


*N 

2a 

o 


r "n 

0 

B/2 (A-l) B/2 


a ic 


1 

B/2 (A-4) B/2 


a 2c 


0 

B/2 (A-9) 


a 0 

3c 


0 

(A-l) B/2 

4 

a ls 

► =4 

0 ' 

B/2 (A-4) B/2 


a 2s 


0 

B/2 (A-9) B/2 


a_ 

3s 


0 

B/2 (A-16) 


a, 

L 4s J 


0 


In this particular case a^, a 2g , a , a^ g =0. But of more 

importance the harmonic sines and cosines may be solved for separately. This 
is a consequence of the absence of velocity dependent terms and does not 
occur in general. 

A property that does generally occur for the helicopter equations, how- 
ever, is the fact that the forcing functions occur only in the first harmonic. 
In the case at hand the forcing functions exist only in the first two rows, 
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i.e., for the mean value and first harmonic. The equations of the higher 
frequency components are not forced; so high frequency coefficients may be 
solved for algebraically in terms of low frequency coefficients > so that 
ultimately the high frequency coefficients may be expressed solely in terms 
of the two forced coefficients. Choosing the set of equations to the third 
harmonic allows the property to be illustrated^ as follows: 


A/2 B/2 


2a } 
o 


^ 0 \ 

B/2 (A-l) B/2 


a icl 

— 

1 1 

B/2 (A-4) B/2 


a 2c 

- 

0 I 

B/2 (A-9)_ 


a_ 

3c j 


1 0 / 


Eliminate a„ and a^ from the equation as follows: 

3c 2c 

[ 0 /2] a 2c + [A-9] a 3c % 0 

if a 4c ~ 0 ( if a 4 c is neglegible) 


a 3c = - [ A -9]' 1 [B/2] a 2 C 

and 

[B/2] a lc + [ A -4 ] a 2c + [ B/2 ] 

[B/2] a lc + [[A-4] - [B/2] [A-9]' 1 [ B/2 ]] a^ 
a 2c = ” [l> 4 ] - t®/ 2 ] [ A_ 9 ] _1 [B/2]] _1 [B/2] a lc 


The equations may then be reduced to two unknowns: 


A/2 B/2 

B/2 [[A-l] - [B/2][[A-U] - [B/2] [A-9]' 1 [b/2]]" 1 [b/2]] 


2a 


lc 


0 

1 


It is interesting to inspect the detailed structure of the (2,2) element 
[A-l] - [B/2] [[A-4] - [B/2] [A-9]' 1 [B/2]]" 1 [b/ 2 ]1 
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It may be expanded to any number of harmonics "n" by analytic 
continuation and takes up the general form: 

[[A-l] - [b/2][[A-4] - 0/2] [[A-9] 

... - [b/2][[a - (n-1) 2 ] - [B/2][A-n 2 cutoff J- 1 [B/2]]- 1 [b/ 2 ] ... 

] ^ [B/2]] _1 [B/2]] 

The value of the (2,2) element converges very rapidly as additional 
harmonics are included. Generally speaking three or four harmonics give 
adequate accuracy to the mean and first harmonic term and sometimes to the 
second. As an example, solve for the harmonic coefficients of 

w + (A + B cos t) w = cos t 

for A = 2.0, B = 1.0 

First evaluate the (2,2) element with two, three, and then four har- 

monics at cutoff. 


Number of Harmonics Value of 

Considered (2,2) Element 


2 


1.125 


3 1.1272 

4 1.12727 


The 

equations to 

solve become 


A/2 

B/2 


2a 

0 


B/2 

(2,2) 


a_ 

lc 

and 


a 

0 

= 

.285 


a ic 

= 1 

.14 


2 by 2 matrix system as follows: 

0 I 

1 i 


With the mean and first harmonic known, the expanded set of equations 
may be employed to determine the next few higher harmonics approaching the 
n cutoff limit, but of course not exceeding it. 
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.2906 


= .022 8 

3c 

The solution of the forced oscillation of the equation then becomes 

w = - .285 + l.l4 cos t + .29 cos 2t + -023 cos 3t 

and is shown in Figure 13. 

It is interesting to note that even though no mean force exists a mean 
displacement of the response occurs. Physically speaking, this is a con- 
sequence of the negative excursion of the sinusoidal force acting on a reduced 
value of stiffness. 


This process is important in explaining the large effects on the mean 
aeroelastic derivatives of including the harmonic coefficients in the ver- 
tical motions differential equations at advance ratios greater than unity. 


Why there is no subharmonic response. - In postulating a form for the 
forced response of the system with periodic coefficients the foregoing 
discussion has shown how the higher harmonic coefficients of the response 
are coupled together so that, in fact, an infinite series of harmonics exists. 


It is logical to consider whether the system does respond also in a 
subharmonic fashion. This question is answered in a straightforward way. 
Assume the series to contain subharmonic terms in addition to the superhar- 
monic so that it appears as follows: 

w = .... ai sin it + ai cos it + a + a n sin t - 1 - a, cos t + 

p:S ? ? o Is 1c 


w 

w 


. . . . p a j t/2 
• • • ""yt sin t/2 


- i a /_ sin t/2 + 0 + a., cos t 
^ c/2 7 Is 

o , cos t/2 + 0 - a n sin t - 

!■ a c / 2 ' Is 


- a_ sin t +. . 
lc 


a n cos t + . . 
lc 


Substituting displacements, velocities, and accelerations in the dif- 
ferential equations adds the subharmonics to the existing terms. They 
become : 
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Figure 13 . Steady Oscillatory Response of a Forced 
System With Periodic Stiffness 







. | a S/ /2 Sin " 17 a c ^2 cos V 2 • • • • 

....A Si s j2 sin V 2 + A a c /2 C ° S 


.... B a s ^2 cos ^ sin V 2 + B a c/2 COS ^ COS ^/ 2 ••• 


and expanded they are as follows: 
1 


a s/2 Sin t / 2 “ 1 a c /2 COS V 2 
+ A a g/ /2 S ^ n V 2 + A a r*/P COS V 2 


c/2 


.... + B a s y 2 sin !*5t - ^ sin t/2) 
+ B a c ^2 (2 cos l*5t + Jr cos t/2) 


With 

all terms assembled together, 

the matrix becomes 

: 



A/2 

0 

B/2 0 

| 



(2a " 

O 


r o" 

0 

(A-l) 

0 B/2 




a -, 

Is 


0 

B/2 

0 

(A-l) 0 

1 

' I 



a ic 


1 

0 

B/2 

0 (A-4 ) . 

1 



a 2s 


0 

• 

• 

• . . 

1 


1 

• 

, = < 

0 

► 

• 

_ _J_ 

. * • 




* 


0 




Hi’ 

0 ... 


a , 
s/2 


0 




‘ 0 

1 

( A -ihJr) 


a / 
c/2 


0 




1 

• • 




0 

y. J 


It is evident that the subharmonics are not coupled to the superharmonics 
in any way. They could only be excited by being forced separately, and in 
this problem that subharmonic forcing is absent. 

It may be concluded that a linear system forced at the period of its 
differential equation periodic coefficients possesses superharmonic response 
by no subharmonic. 
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It should be noted that the real system may not be strictly linear, and 
in that case, subharmonic response is a distinct possibility. 
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Three- and Four-Blade Rotors 


The basic method developed in the general discussion will now be applied 
to rotors having three and four blades. The solutions will be similar ex- 
cept for the inclusion of the rotor reactionless mode in the four-blade 
solution. This mode, of course, does not exist for a three-blade rotor. The 
method that follows will be developed primarily with a four-blade rotor in 
mind, but may be adapted to the three-blade rotor simply by dropping the 
matrix rows and columns associated with the reactionless mode and converting 
the 4P terms to 3P* 

The nine degree-of-freedom equations in stationary coordinates are as follows 

Wf's'l +[[■>] + [a$9} + [[S] + [Ap]]|B[ =[fCFl + [A F l]{ nF | + | M ( (l) 

where the matrices contain the following terms: 
inertia 

mechanical damping and gyroscopic 
aerodynamic stiffness derivatives 
structural stiffness 
aerodynamic damping derivatives 

contains a centrifugal collective flapping force term due 
to rotor precone 

aerodynamic derivatives due to fixed angles and blade geometry 

forces and moments applied to rotor and gyro 

In addition to mean values, the aerodynamic derivatives contain terms 

which are functions of bfit, 2b fit, NtOt, where b equals the number 

of blades. The derivatives that are functions of 2tflt, , Nbftt will be 

neglected, since their values are small compared with those of the mean and 
the b per rev harmonic terms. 


[i] 

O] 

[A P ] 

[s] 

[ A (d 

[cf] 

[M 

M 
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The four-bladed rotor derivatives also have 2P, 6P, . .., (2 + 4 n)P 
harmonic terms. Those higher than 6P will be neglected as they are rela- 
tively small. 

Expanding the aerodynamic derivatives in terms of their harmonics while 
combining structural and mean aerodynamic terms, equation (l) becomes: 

[i]H1e o Me 2c ]cos 20t + [E^Jsin 2fit + [E^jcos + [E^J sin 

+ [ E 6 c h OS 6nt + [ E 6 s b in ^t]{e}+[[F o ] + [F 2c ]cos 20t + [F 2s ] sin 20 1 

+ [ f 4 c ] cos + [ F u s h in + [ F 6 c b os 604 + [ F 6 s b in M = [[ G 0 ] 

+ [G 2c ]cos 20 1 + [G 2s ]sin 20t + [G^Jcos IjOt + [G^Jsin ¥5t + [G^ c ]cos 60 1 

+ [ G 6 s ] sin + | M ( ( 2 ) 

(For three blades, equation (2) would contain only mean and 3P terms.) 

Equation (2) must be solved for 8, which is a function of time. Let 8 
be represented by a Fourier series: 

3 = 0 + 8 n cos 20 1 + 8 0 sin 2fit + 6) cos Wt + 8* sin Wit 

o 2c 2s 4c 4s 

+ 8g c cos 60t + 8^^ sin 60 1 + 8g c cos 80 1 + 8g s sin 80 1 

+ cos ICOt + 8_, _ sin 10 0t + ... 

10c 10s 

The series will be truncated at the l4P term. Coefficients of the higher 
harmonics will be negligible. 

Equation (2) will be solved according to whether the swashplate is fixed 
or free and the vehicle is free flying or has its shaft fixed. 

Fixed-shaft, fixed- swashplate solution. - The vector of displacements 
is as follows : 
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where is the differential collective mode (also known as the reactionless 
or the self-balancing mode). 

For the fixed-shaft f ixed-swashplate solution: 

0 = 0 = 0=0 

¥ = $ = $ _ o 


z = z - z - 0 

e = e = o 

^ =0=0 


6 and 0 may take steady nonoscillating specified values. The body de- 
grees of freedom are locked out of the equations. Thus the solution becomes 
applicable to a vehicle rigidly mounted in a wind tunnel and controlled by 
displacement inputs to the swashplate. 



may be divided into upper and lower portions: 
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where the given value of < (3 u V is: 



It is required to solve for the displacement vector { }• • However 

{ } is complicated by the reactionless mode. Its coefficients in the dif- 
ferential equations include 2P and 6P components as well as the 0 and 4p 
expected for the four-blade rotor. But since it is to be a steady oscillatory 
response, each blade in turn must execute identically the same motion as the 
blade ahead of it. It is clear, therefore, that hub moments, thrust, swash- 
plate moments, -in fact all data measured relative to the stationary axes-must 
contain only 0, 4-P, -8P, etc, harmonics. This means that the reactionless mode 
must exercise vigorous 2P, 6P, etc, motion in rotating axes to attain a 
stationary blade pattern in stationary axes. It should be possible, therefore, 
in stationary axes, to represent the modal motions selectively as follows 
rather than to employ all harmonics of all terms simultaneously. Let 
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where 


{ B *uf = l 6 ul = { 6 u \ + IW COS + \ 6 uUs I Sin ^ 

o 

+ I 6 u8c^ COS 311 + I ^u8s I sin 801 + 

and 

{ 6 dl = { 6 d2cl COS ** + { 6 d2sl Sin ** + I 6 d6cl cos 
+ Pd6sl Sin 601 + *••• 

By similar reasoning, the last column and bottom rows of the [Ap] and 
[A P ] matrices of equation (l) will contain only 2P and 6P terms of significance 
(except for the on-diagonal term, which will contain only mean, 4p, and 8P 
terms of significance). The bottom row of M in equation (l) will also 
contain only 2P and 6P terms of importance. 

Recalling that 



equation (2) may now be rewritten in partitioned form: 
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The shaded areas contain 2P and 6P and 10P, .... only, all the rest 
contain 0, Up, 8P, .... 


This is permissible because I = I = I = I = 0 and the products 

1 J Jjl 2) 

of (2P, 6 p ) x ( 2P, 6P ) = 0, 1+P, 8 p, 12P 

(2P, 6p ....) x (0, UP ....) = 2P, 6p, 10P .... For the present solution, 
the equations simplify to 
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Recall that the vector of applied moments in equation (4) is: 



| M } may be redefined to represent reaction loads (aeroelastic response 
loads). Thus the vector may be moved to the left-hand side of equation (4) 
with no change in sign. Additionally, by moving the terms involving the 
input vector | 8 U | to the right-hand side, equation (4) may be rewritten: 



Equation ( 5 ) may be separated into uncoupled upper and lower halves. 


The lower half may be written: 



This equation is the basis for determining the displacement vectors. 
Once they are known, the aeroelastic response loads can be found from the 
upper half of equation (5), which is as follows: 

| M f = " “ F E 12 E 13l f 6 U } " [ F 12 E 13~J I"! 

L J { dj L J l dj 

+ C G i]{M " ["iiJIM 

where the complete right-hand side of the equation is known. 
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Solution of forced equations of motion: - The problem now is to solve 
for 6 U and 6^, starting with equation (6). 

= ^ u ) + 6 u cos 4flt + 6 sin bClt + 6 cos 8^t + 6 sin 8Ht + 

° 4 c 4 s 8c 8s 

6 = 6, cos 20 1 + 6 sin 20t + 6 cos 60 t + 6 sin 6ft t + . .. 

d d 2c d 2s d 6c d 6s 


The procedure is tedious but straightforward to solve for the values of 
the coefficients of the series. Each matrix is expanded in its respective 
harmonic parts. The series representing the solution is also expanded and 
products of matrices and solutions are found and all terms of common frequency 
and phase combined and a truncated solution found. This is shown in detail 
in the following paragraphs. 


First, abbreviate the notation for clarity: 
Let 

20 1 = t 1 


and 


2 O = Q J 


6 u = B 


subcripts 

C 1 = 2 c 
s"*" = 2s 


6 =6 
d 


Therefore, 


6=0 + 0 O 1 cos 2t"^ + B_ 1 sin 2t^ + 6,1 cos 4t^ + 6,1 sin 4t^ + 

o 2c 2s 4c 4s 

6 = 61 cos t^ + 6 1 sin t^ + 6 1 cos 3t^ + 6 1 sin 3t^ +, 

lc Is 3c 3s ^ 


and let the matrices be abbreviated. 


E 22 

= M 

F 22 

= R 

E 23 

= N 

F 23 

= S 

E32 

= P 

F 

32 

= T 

E 33 

= Q 

F 

33 

= U 


[ G 2 : - F 2l] 

[ G 3 ! - F 3 i] 


V 


W 
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Then equation (6) becomes: 


p -1 

/ .. \ 

|- 

/ • \ r 

*22 ° 

1 8 1 

M N 

i e 



J j + 


t 


° I 33 

6 

P Q 

1 * 


_ _ 

' * 


' L 



The two components of the vectors and their derivatives are as follows: 

8 = 8 q + B^l* cos 2t 1 + s ^ n ^t 1 + B| +c 1 cos 4t^ + 1 s ^ n ••• 

B = -2D 1 B 1 sin 2t 1 + 2ft 1 0 1 cos 2t^ - 4fl 1 8| 1 sin 4t^ + 4ft 1 0, 1 sin 4t^ 
2c 2s 4c 4s 

2 2 4 

B = -4ft 1 8 2c 1 cos 2b 1 - 4ft 1 ®2s 1 sin ^4 C 1 cos ••• 

p 2 

l6n 8, 1 sin 4t 
4s 


6 = V COS ^ + ^ls 1 sin 11 + ^c 1 cos ^ + 6 3s 1 sin ^ ••• 

6 = ” q16 - 1 1 sin t 1 + n 1 ^ 1 cos t 1 - 1 sin St 1 + 1 cos St 1 

-Lc is 3c 3 s 

2 2 2 2 
6 '= -0 1 6^ cos t 1 Q 1 6 lg i sin t 1 - SO 1 6^1 cos St 1 - gfi 1 6 1 sin St 1 

and the matrices become 

M = M + M 1 cos 2d 1 + M 1 sin 2t^ 
o 2c 2s 

N = N lc l cos t 1 + N lg l sin t 1 + N^l cos St 1 + N^l sin St 1 

p = P lc ,l cos t 1 + P lg l sin t 1 + P^l cos St 1 + P sin St 1 

Q = Q q + Qp c ^ cos 2p + Q 1 sin 2il 

For the R, U equations,, substitute R, U for M, Q 
For the S, T equations, substitute S, T for N, P 


The work becomes two matrix equations in the two unknown vectors, 8, 6 
when expanded, in their harmonic components. 
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Treating the 3-Force equation first, matrix by matrix, yields the 
following expansion from equation (8): 

2 2 2 
I 22 B 2c 1 COS 2tl * I 22 B 2 s 1 Sin 2t± ” V ^ I 22 B 4c 1 G ° S ^ 


t M + M„ 1 cos 2z~ 
o 2c 


2 2 
■ 4CJ 1 i 2 2 B 2s 1 Sin 2t± ” 1601 I 22 B 4c 1 COS 1+t ' L 

l 2 1 

160 I 22 B 4s 1 Sin ^ 

T f 


M„ 1 sin 211 -2fd|3 0 1 sin 2 d + 2fds o 1 cos 2 d 

2s J 2c 2s 


+ I N, 1 cos d 
lc 


J I 

-l40 1 B 1+c l sin bt 1 + 40 ^l cos 1+tl 

4- N 1 sin t^" + N 0 1 cos 3t"^ + N 0 1 

Is 3c 3s sin 


)s t^* + N-, 1 sin t^" + N 0 1 cos 3t"^ + N 0 1 . 0 . 1 r>lc 

Is 3c 3s sin 3t -wo 

\ 

> 1 cos t 1 - 3n 1 6„ > 1 sin 3-t 1 + 3fdfi 0 1 cos St 1 

Is 3c 3s 

4 1 cos 2t^ 4- R 1 sin 2t | |8 + B_ 1 cos 2t^ - 

>/-> O e 1 I I ^ On 


+ 0^6 l cos t"^ - 3^ 1 6 0 1 sin 3t^ + 30^6 
Is 3c 

+ R_ 1 cos 2t 1 *f R 1 sin 2t x 8 + 8 0 

d 2c 2s o 2 


> 1 sin t 

lc 


^ 4- R 1 sin 2t' L 
2s 

4- B^ 1 cos bt 1 - 


+ 8^ 1 cos 2t^ 
o 2c 


4- 1 sin 4t^" 


4- S_ 1 cos t^ 
lc 


cos t 1 4- S n 1 sin t 1 + S„ 1 cos 3t^ 
Is 3c 

5 1 sin d + 6 1 cos St 1 + 6 1 sir 

Is 3c 3s 


= V + V 0 1 cos 2t 1 
o 2c 


>, 1 sin d + S_, 1 
Is 3c 

+ 6 l cos 3d 

3c 


V 0 1 sin 2t 1 
2s 


S 1 sin 3 d fi. 


• 04-1 4- 

> 1 sm 3t + 

3s 


+ b 2s i s ^ n ^d 


X ] [he 1 cos ^ 



Treating the 6-Force equations yields, by expanding equation (8): 

2 2 2 2 

I ( Q 1 6 1 cos d - fd 6 1 sin d - cfi' 1 ' 6 l cos 3 d - $fd 6 1 sin 3 d) 

33 lc Is 3c 3s 

T 1 0 \ / T 


^P lc l COS d + I 

^5 1 B 0 1 COS 2t 1 
2s 


„ , . ,1 
P, 1 sin t 
Is 

1 - *«\ 1 1 


re 1 cos 2 d - 4n 1 e 4c i sin 4d + 1 

Q + CL 1 cos 2 d + Q_ 1 sin 2 d) (-d~6 1 
o h2c 2s lc 


sin 2 d 


.s , - 1 , 


Ln d + P_ 1 cos St 1 + P_ 1 sin 3 d) /-2fdp sin 2 d 
3c 3s ' \ 2c 

j^l sin 4t 1 + 4fdp^ 1 cos 4t 1 ) 
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( 


T 1 cos t 1 + T.. 1 sin t 1 + T 0 1 cos 3t'*’ + T 1 si 
lc Is 3c 3s 


B o + 


in St 1 ) 

+ 13^1 cos 4-t 1 + Bj^l sin 4-t 1 ....) 

+ U 2c l cos 2t 1 + U 2s l sin 2t 1 ^ cos t 1 + 6.^1 sin t 1 

+ 6 l cos 3t 1 + 6 l sin St 1 ...) 

3c 3s / 

L W 1 cos t 1 + W n 1 sin t 1 + W_ 1 cos St 1 
lc Is 3c 


fi 1 cos 2-t 1 + 1 sin 2t 1 

2c 2s 


+ W 0 1 sin 3t 
3s 


f ft ^ 


Z 

e 

l 0 


(io) 


Expanding the terms of the 8- force and 6 -force equations and employing 
trigonometric identities to reduce all terms to harmonic sines and cosines 
of the first degree allows a set of simultaneous equations relating the 
coefficients of the Fourier series to he assembled. Each equation of the 
set contains terms of one harmonic component only. Dividing each equation 
by its harmonic yields a set which can be truncated and solved for the 
coefficients , algebraically . 

The set of equations may be written in the following order. In all 
following work the Fourier series is truncated at l4P. This allows the vec- 
tors to be described in terms of the following harmonics: 


B: 

B , 

o’ 

6: 

6 

lc 


6 

Is 


2s 

3c 


6 

3s 


6 

5c 


0 6c’ B 6 s 

6 , 6 , 
5s’ 7c’ 


6 

7s 
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Combining some of the A_ matrices yields: 


r~ i 
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C, X ^ 


r 

A 1 !-! 2 - 


B i4 

0 

0 

0 

0 


(K) 0 


(F u4 

B 21 1 B 22 

I 

B 23 

B 24 

B 25 

0 

0 

0 


<44 


<44 

B 3i 1 B 32 

B 33 

B 34 

B 35 

B 36 

0 

0 


(« u ) 2 


(Fu ^2 

B 4l ] \2 

B 43 

B 44 

B 45 

B 46 

B 47 

0 

•i 

(6 ) 

V 1/3 

► = -< 

(f ) 

V L ; 3 

0 l B 5 2 

B 53 

B 54 

B 55 

B 56 

B 57 

B 58 




0 

1 

0 I 0 

1 

B 63 

B 64 

B o5 

B 66 

B 67 

B 68 


<44 


0 

0 1 0 

1 

0 

B 74 

B 75 

B 76 

B 77 

B 78 


< 6 u >6 


0 

0 I 0 

1 

0 

0 

B 85 

B 86 

B 87 

B 88 


(* 1)7 


0 
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M 
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< 4 >_ 


= 

rA 22 

*23 


'"■32 

*33 _ 
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1 

1 — 1 
C\J 
< 

1 



, A 3i j 



ho 1 



6 1 
Is 

► 


L - 



etc . 


etc . 


(«u). 


^c 1 


^s 1 


etc . 
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By virtue of the zero forcing functions in the 4, 5, equations, 

the matrix may be reduced even further for calculation of the harmonic re- 
sponse, by eliminating the ( 6 U )^ 5 ( 6^)5 ? (^u)g ? terms in the 

same fashion as shown for the one degree of freedom example. The reduced 
matrix then becomes l 


l4l 




A n 

B 12 

B 13 

B 14 ‘ 


'< 8 G 0 ' 







B 21 

BC 22 

BC 23 

BC 24 


(‘it 



< F L>1 




B 31 

BD 32 

BD 33 

BD 34 


(« u >2 

► 

= M 

1 

(F u ) 2 

► 

(13) 


B 41 

to 

BE 42 

BE 4 3 

BE 44_ 


(6 ) 

w L 3j 



'J r A 



and relationships of 

the following type allow the calculations 

of the higher 

harmonic components : 










( 6 u) 4 

= 

C 52 

(6 ) 
v i/l 

+ 

C 53 ^ u ^2 

+ 

C 54 

(6 ) 
v 1/3 



<V 5 

= 

C 63 

(«u) 2 

+ 

C 64 (6 l } 3 

+ 

C 65 

( 6 u) 4 

(14) 


< 6 »>6 

= 

C 74 

(6 ) 

V 1/3 

+ 

C 75 (6u) 4 

+ 

C 76 

(6 ) 

V 1/5 


(6 ) 

' 1 7 

= 

C 85 

( 6 u) 4 
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C 86 (6 : 


+ 

C 87 

(6u) 6 


where the ! 

B C, BD 

, BE 

& C. . 

10 

are determined 

from 

values 

of A . . . 
ij 



With the responses of the degrees of freedom found, it is possible then 
to find the thrust, hub moment and swashplate mean and oscillating forces, 
blade loads, and blade tip motions. 


Hub moments, thrust and swashplate moments: - With the motions of the 
degrees of freedom known it is possible to calculate the forces and their 
harmonic components |m£ from equation (7) 






[°i] IM 


(7) 


and these equations may be expanded in terms of all harmonic components and 
as before these may be arranged in order of the harmonics. Grouping ac- 
cording to harmonics the equation takes on the following form: 
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The form of the D 1 s is similar to the A's of the A. . matrix. 

Tip path motions. - With the shaft and swashplate both fixed against 
tilting oscillation, the only motions that occur in the rotor system are the 
deflection of the blades. This section, therefore, sets out to show what 
these deflections are over the range of advance ratio and flapping stiffness 
ratio P for which these analyses are expected to be valid. 


The calculations were performed employing the Lock number y = 4.57 of 
the 33-foot rotor. The nondimensionalized blade tip deflection 6/R 3 however, 
should be valid for any rotor, with any number of blades, if it has the same 
value of Lock number. 


Figures l4 through l6 show the variation of tip deflection as the blade 
travels the azimuth for those familes of conditions ^ = .5? 1.1, and 1.7 
with members of each family of flapping stiffness ratios of P = 1.5, 2.0, 
and 5*0. 
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Figure l4 shows the rate of change of tip path deflection with respect 
to Figure 15 the rate with respect to 0 1o , and Figure l6 shows the 


Is' 


residual tip path deflection occurring at 9 n = 0 =0 and caused by the 

lc Is 

blade twist, rotor precone, and the small collective pitch, © =1*5 , 

• 75R 

present in the tests of the 33-foot rotor. 


Mean aeroelastic derivatives. With the fixed-shaft, fixed-swashplate 
equations solved for the motions of the blades, it is possible to calculate 
the shaft-transmitted forces, blade airloads, and flapping bending moments 
and the input data required for the in-plane equations of motions. This 
section of the report is concerned only with the mean values of the aero- 
elastic forces transmitted through the shaft, the mean aeroelastic derivatives. 

Mean aeroelastic derivatives of the rotor are commonly used in arti- 
culated rotor helicopter stability analyses at low advance ratio. For 
gyroscope stabilized hingeless rotors, especially at high advance ratio, it 
is usually not satisfactory to separate overall body motions from those of 
the rotor gyroscope system. This makes the concept of mean aeroelastic 
derivatives of limited usefulness in these applications. Certain of the 
derivatives, however, are easily measured in wind tunnels and provide an 
excellent reference by which to judge analysis methods. It is for this 
reason that they are treated here. 


Theoretical and experimental values of rotor mean aeroelastic deriva- 
tives are shown for the 33-foot three-blade rotor and the 7.5-foot four-blade 
rotor. First the derivatives of hub moment, swashplate moment, and thrust 
with respect to cyclic pitch for the 33-foot rotor, are discussed. 

The theoretical variation of hub moment coefficient derivatives with 
advance ratio u and flap frequency P are shown as vectors in the X, Y plane 
in Figures 17 and l8. The derivatives are applicable to rotors of any 
size, numbers of blades, and blade geometry, provided they have the same 
Lock number (V = 4.57) as the 33-foot rotor. Two sets of curves are shown. 

The heavy lines show the effects of including the harmonic components of the 
aerodynamic coefficients in the differential equations. The light lines 
show the effect of leaving them out, as in the common practice with arbiculated 
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Figure 18. Hub Moment Lateral Cyclic Pitch Aeroelastic Derivatives, V = 4.57 




rotors at low advance ratio. It is interesting to note that at advance ratios 
well below unity, **<1.0, there is little difference between the two sets of 
curves. This seems to justify the common practice of leaving out the 
periodic coefficients. At advance ratio greater than unity, *i>1.0, however, 
the effects of the harmonic components are very large. A large bump de- 
velops in the vicinity of P = 2.0, in the derivatives with respect to 0 , 

as the advance ratio increases above ** = .8 so that the mean derivatives 
become much larger with the periodic coefficients included. With P<1.8 the 
effects are large, but the derivatives with periodic coefficients become 
smaller than their counterparts without periodic parts. 

Comparisons of theoretical mean aeroelastic derivatives with the experi- 
mental values for the 33-foot rotor are shown in Figures 19 through 22. 

In these plots the components of the derivative vectors are displayed 
versus flapping frequency ratio P for tested values of advance ratio. 

Small corrections have been made to the experimental points to bring 
them to common values of advance ratio. There is general quantitative 
agreement between theory and experiment but the fine structure of the 
variation with P exhibited by the experimental data is not mirrored in the 
theoretical results. This may be due to too restricted a mathematical 
representation of the participating vibratory modes in the theory. 

In Figure 21, experimental hub pitch moment aeroelastic derivatives due 
to longitudinal cyclic pitch for the 7*5 foot diameter 4-blade rotor are in- 
cluded with the 33 foot rotor data for comparison. The derivatives are inter- 
polated to produce values at advance ratio ** = 0.5? 0.8, 1.1, and 2.0 and 
are shown at four values of flap frequency ratio. It is interesting to note 
that the 7-5 foot 4-blade rotor data at ** = 0.5 and 0.8 compares very favor- 
ably with values measured on a rotor of more than four times the diameter, 

that possesses only three blades and has about half the solidity. 

Figure 21 also shows a significant effect due to the blade first flap- 

wise radial mode shape. The 7*5 foot rotor configuration with the stiff 
flexure gives greater values of pitch derivatives, for the same **, y, and P, 
than does the configuration with the less stiff flexure. The effects of the 
modal differences seem to increase with advance ratio. 
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FLAP FREQUENCY RATIO, P 


Figure 21, Mean Aeroelastic Derivative of Hub Pitch Moment With 
Respect to Longitudinal Cyclic Pitch, V = 4.57 
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Swashplate moment coefficient mean aeroelastic derivatives with respect 
to cyclic pitch for the 33-foot rotor are shown plotted in vector form in 
Figures 23 and 24. It should be noted that in addition to being 
divided by rotor solidity o 9 as were the hub moment coefficient derivatives , 
the swashplate moment coefficients are divided by the blade sweep angle A 
and factored by the swashplate mechanical advantage k. These further 
factors serve the purpose of making the derivatives independent of A and k 
at low advance ratio. At high advance ratio the derivatives become a weak 
function of sweep ratio, ^ (or rr ) 5 because of reverse velocity 
effects . 

Swashplate mean derivatives variations with P and ^ at Y = 4.57 are 
similar to those of the hub moments. The effect of the cant angle = 60° 
between the blades and the swashplate is apparent in the approximately 60° 
counterclockwise skew of the derivatives relative to the hub moment deriva- 
tives. As with the hub moment derivatives, deletion of the harmonic com- 
ponents of the coefficients in the differential equations has little effect 
at advance ratio less than unity and a large effect at values greater than 
unity. 

Comparisons of theoretical swashplate moment coefficient derivatives, 
calculated with harmonic components included, with values measured 
experimentally on the 33-foot rotor, are shown in Figures 25 through 

28. 

As with the hub moment derivatives the vector components are shown 
versus flap frequency P at particular values of advance ratio Again the 
agreement between theory and experiment is quantitatively good but the fine 
structure of the experimental data is not seen in the theory. 

The variation of the thrust coefficient C^/o with cyclic pitch is 
shown in Figures 29 and 30 over the range of P and |i for which the 
theory is expected to be valid. The theory includes the effects of the 
harmonic components of the aerodynamic derivatives. Agreement between ex- 
periment and theory is better for the longitudinal cyclic pitch because of 
its greater effect. 
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Figure 26. Mean Aeroelastic Swashplate Roll Moment Derivative With 
Respect to Lateral Cyclic Pitch, V = 4.57 
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Figure 30. Mean Thrust Aeroelastic Derivative With Respect to 
Longitudinal Cyclic Pitch. V = 4.57 


The residual shaft force coefficients for the 33-foot rotor, or those 

existing at 6, = 6. = a = 0 and 0 = 1.5°, are needed to reproduce the 

^ lc Is . 75R 

mean states in which the rotor was tested. Mean shaft force states may be 
obtained by combining residual force with the forces produced by the two 
cyclic pitch components of the test conditions desired. 

The residual mean force coefficients, hub moment, swashplate moment, 
and thrust were caused by blade twist, rotor precone, and the small col- 
lective pitch angle 0 ^ = 1.5°. The theoretical variation of them with 

flap frequency P at the specific values of M* tested are shown in Figures 
31 through 35 • Experimental values are also shown. Good agreement was not 
expected and did not occur. The reasons for the poor agreement are thought 
to be the following: 

a. Flow distortion caused by the body 

b. Centrifugal flattening of precone not adequately described by one 
parabolic mode 

c. Blade bending in the second flap mode due to blade twist aerodynamics 

d. Induced inflow due to blade twist aerodynamics not accounted for 

It should be noted that some of the residual forces could not be ob- 
tained experimentally without overloading the rotor. They have been obtained 
by extrapolating best fit plane (or rms plane) data to the zero cyclic con- 
dition. To avoid giving a misleading impression, therefore, the residual 
forces in conjunction with the mean aeroelastic derivatives have been em- 
ployed to reproduce a mean state of interest - the condition of hub moment 
trim (or zero hub moment). Cyclic pitch angles needed to trim the 33-foot 
rotor hub moments to zero are shown in Figures 36. The actual test values 

of and 0, were centered about the trim values and generally did not 
lc Q Is 

exceed 3 or 4 amplitude deviation from them. 

Residual harmonic aeroelastic forces are discussed later and comparable 
cyclic pitch values may be found to trim or minimize vibratory loads. The 
7.5-foot four-blade rotor hub moment and thrust coefficient aeroelastic 
derivative theoretical values, divided by solidity, are very similar to those 
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FLAP FREQUENCY RATIO, P 


Figure 31, Residual Mean Hub Pitch Moment, Y = 4.57, 

B o = 2.25°, e t R = -9.43°, or = o, e >75R = 1 . 5 ° 
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Figure 35. Residual Mean Thrust Force Due to, Y = 4.57, 0 = 2.25°, 

e t R = -9.43°, c* = 0, and 8^ = 1.5° ° 
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of the 33-foot three-blade rotor. The difference, about 10 % 9 is due to the 
larger Lock number, V = 5.0, of the four-blade rotor compared to that of the 
three -blade rotor, Y = 4. 57. 

For this reason the theoretical derivatives hub moment vector diagrams 
with respect to and 0^ are not repeated. In this section the hub moment 
derivatives components of the vector diagram are plotted versus advance ratio 
for the test values of flap frequency P. 


The following derivatives are plotted: 
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and comparison with experimental values is shown in Figures 37 through 44. 

It should be noted that residual forces for the cantilevered blades with 
no precone are theoretically zero. 

The 7-5 foot 4-blade rotor test procedure and experimental data are 
discussed in detail in Reference 17* 
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Oscillatory aeroelastic derivatives . This section of the report is 
devoted to the presentation of calculated and experimental values of shaft 
and swashplate - transmitted oscillatory forces for the fixed-shaft, fixed- 
swashplate condition. Shaft shear force is a consequence of in-plane motion 
and is discussed in the section on solutions of in-plane motion equations. 

In the section on mean aeroelastic derivatives it was possible to 
discuss three- and four-blade rotor results together, when suitably nondi- 
mensioned, since the results did not depend on the number of blades in the 
rotor. For oscillatory forces transmitted through the shaft and swashplate 
however, the number of blades determine the frequency and amplitude of the 
result. For this reason three-blade and four-blade rotors are treated 
separately. 

The theoretical and experimental results for the 33-foot three-blade 
rotor are presented first. These are followed by the results for the 
7.5-foot four-blade rotor. 


The theoretical harmonic response aeroelastic derivatives of hub and 
swashplate moment coefficient with respect to 0^ and 0 , divided by solidity, 
are presented in vector diagrams for the 33-foot three-blade rotor. In station- 
ary axes the hub moment oscillations contain only 3P oscillations but these are 
conveniently thought of as having been produced by 2P advancing and 4p regres- 
sive oscillations relative to coordinates rotating with the rotor. The 2P ad- 
vancing contributions to the 3P stationary axis hub moment derivatives are shown 
in Figures 45 and 46, and over the ranges of P and ^ for which these analyses are 
expected to be valid. Figure 46 showing the variations of the hub moment coeffi- 
cient vectors, divided by solidity, with repect to 9 shows that at low values 

-L S 

of advance ratio, say u = 0.5 for example, and at low values of rotor blade 
stiffness ratio, say up to P - 1,4, that the 2P contribution to vibratory force 


is small compared to the mean hub moment produced, as seen in Figure 17. It 
... , i. , ^ , I r'Cm/ cr*l 


consists at P = 1.4 of nose-down pitch moment of about 


* = -.05 


whereas the mean hub moments generated at the same condition are dCm/o 


.3^ or not much more than about -j^-of the mean hub momen?. 
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Figure 45. Oscillatory Hub Moment 2P Aeroelastic Derivative With Respect to Lateral Cyclic Pitch, 
V = 4.57, Vectors Shown at ^ = 0 
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Figure 46 . Oscillatory Hub Moment 2P Aeroelastic Derivative With Respect to Longitudinal 
Cyclic Pitch, y = 4,57 Vectors Shown at if = 0 



At P values typically employed in conventional hingeless ratio helicopters , 
i.e., P = 1.2, the vibratory contribution would be negligible. 


It is interesting to note, still at M .5? that if P is allowed to 
reach the neighborhood of 2.0 so that the first flap mode is in resonance, 
the vibratory contribution can become larger than the mean hub moments 
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Should the blades be made very stiff, say P = 10, then the vibratory 
component reduces somewhat but not to as low a value as occurs at low P 
values . 


The same trends continue as advance ratio, p, increases except that the 
vibratory contribution grows more rapidly than the mean moment capability. 

The 2P vibratory response with respect to 0^ c (Figure 45) follows the 
same general trends as did the component with respect to 0 . 

The theoretical hub moment 2P vibratory derivatives with respect to 
0 and 0 , shown in Figures 45 and 46, are replotted separately in Figures 

_LC XS 

47 through 50 versus flap frequency ratio P for the several values of advance 
ratio m* employed in the wind tunnel tests of the 33-f°°t rotor. In these 
plots, comparisons are made with experimental data and it is seen that the 
general trends are predicted well. 

The 4P contributions to the 3P stationary axis hub moment derivatives 

with respect to 9_ and 0_ are much smaller and less well defined than the 
lc Is 

2P contributions . They are displayed in Figures 51 and 52 in general areas 
for each value of advance ratio, to the same scale as the 2P contributions 
and are disassembled and displayed versus P in Figures 53 through 56 for 
the tested values of advance ratio. 

The experimental values are very small at low P values and larger at 
high P values as are the theoretically predicted values, but agreement is 
erratic . 


The swashplate harmonic moment derivatives with respect to 0^ and 
in the 2P contribution to 3P stationary axis moments, display similar behavior 
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FLAP FREQUENCY RATIO, P 


Figure 47. Oscillatory Hub Pitch Moment 2P Aeroelastic Derivative With 
Respect to Lateral Cyclic Pitch, Vector Component at + = 0, 

Y =4.5 
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Figure 48 . 
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Oscillatory Hub Roll Moment 2P Aeroelastic Derivative With Re- 
spect to Lateral Cyclic Pitch, Vector Component at t = o, v = 4.57 
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Figure 50. Oscillatory Hub Roll Moment 2P Aeroelastic Derivative With Respect 
to Longitudinal Cyclic Pitch, Vector Component at ♦ = 0, V = 4.57 
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Figure 51 • Oscillatory Hub Moment UP Aeroelastic Derivative with 
Respect to Lateral Cyclic Pitch, Y = U. 57 . 
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Figure 52. Oscillatory Hub Moment UP Aeroelastic Derivative with 
Respect to Longitudinal Cyclic Pitch* V = U.57- 
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to the hub moment derivatives. As with the mean aeroelastic derivatives, 
the 60° skewing of the swashplate derivative ahead of^the hub values due to 
the cant angle = 60° is evident again in Figures 57 and 58. 

Figures 59 through 62 show the 2P components of the swashplate 
derivative vectors plotted versus P at the tested values of U and compares 
them with experimental data. Although there are not many experimental points 
the theory agrees in sign and order of magnitude with the experimental data. 

The 4P contribution to the 3P stationary axis swashplate derivatives 
with respect to 9^ and 0 are shown in Figures 63 and 64 , The detailed 
variation with P is complicated but within relatively small areas at each 
value of U" The areas are shown. The detailed variation with P and u of 
each 4p component is shown in Figures 65 through 68 compared to ex- 
periment. Agreement is comparable with the 2P swashplate moments. 

Rotor thrust oscillation due to cyclic pitch application is next shown 
theoretically for the 33“P°ot 3-blade rotor - no experimental data is 
available. Thrust oscillations occur in stationary axes as harmonics of the 
product of number of blades and shaft rotation rate tft. Harmonics higher 
than the first are believed to be very small for this rotor; therefore only 
the first harmonic 3P components are shown. 

The sine and cosine components of the thrust coefficient divided by 
solidity, C-j/ct, aeroelastic derivatives with respect to 8 and 9 , are 

-LC _LS 

shown in Figures 69 through 72. It should be noted that at t = 0 the 
position of number one blade is aft. 

It should be noted that the sine component due to 0 and the cosine 
component due to 0^ c both peak at a flap frequency ratio of P = 3.0. 

The residual harmonic forces remaining when 9 n = 0- = & = 0 and 0 

lc Is .75R 

= 1,5 a re shown in Figures 73 through 82. 

The two components, sine and cosine, of the thrust coefficients, divided 
by solidity, C T /a produced by 0 = 1.5* 9.R = -9*^3°* and 3 = 2.25°, are 

shown in Figures 73 and 7b versus flap frequency ratio P for values of 
advance ratio 4 tested. No comparison with test data, however, is available. 
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Figure 57. Oscillatory Swashplate Moment 2P Aeroelastic Derivatives With 
Respect to Lateral Cyclic Pitch. Vectors Shown at 'll = 0, V = 
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Figure 60. Oscillatory Swashplate Roll Moment 2P Aeroelastic Derivative 
With Respect to Lateral Cyclic Pitch. Vector Component at 
* = 0, Y = 4.57 
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Figure 6l. Oscillatory Swashplate Pitch Moment 2P Aeroelastic Derivative 
With Respect to Longitudinal Cyclic Pitch. Vector Component at 
♦ = 0, Y = 4.5 
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Figure 62. Oscillatory Swashplate Roll Moment 2P Aeroelastic Derivative 

With Respect to Longitudinal Cyclic Pitch. Vector Component at 
♦ = 0, V = 4.57 
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63 . Oscillatory Swashplate Moment Up Aeroelastic Derivative with 
Respect to Lateral Cyclic Pitch, V = h . 57 



64 . Oscillatory Swashplate Moment hP Aeroelastic Derivatives with 
Respect to Longitudinal Cyclic Pitch, Y = U .57 
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Figure 65- Oscillatory Swashplate Pitch Moment 4P Aeroelastic Derivative 

With Respect to Lateral Cyclic Pitch. Vector Component at t = 0, 

Y = 4.57 
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Figure 67. Oscillatory Swashplate Pitch Moment 4P Aeroelastic Derivative 
With Respect to Longitudinal Cyclic Pitch. Vector Component 
Shown at i = 0, V =4.57 
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Aeroelastic Sine Derivative With Re- 
ic Pitch, 0 = 4.57 





Figure 71 • Oscillatory Thrust 3P Aeroelastic Cosine Derivative With 
Respect to Longitudinal Cyclic Pitch, y = 4.57 
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Theoretical hub and swashplate coefficient residual 2P and 4p moments 
are similarly displayed in Figures 73 through 82, and compared with ex- 
perimental residual data. 

It should be remarked, in concluding the three-blade case, that the 
residual mean and harmonic rotor shaft forces may be combined with shaft 
force derivatives factored by the appropriate cyclic pitch angles to fully 
reproduce the complete dynamic state experienced by the 33-foot 3-blade 
rotor under any tested condition. 


In summary, a complete rotor shaft and swashplate moment and force 
steady oscillatory state may be assembled for any combination of cyclic pitch 
components, within the linear range, for the 33-foot 3-blade rotor from 
plotted coefficients as follows: 
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Follow the same rules as 

y Cm and C l 
a a 


The nondimensional derivatives are applicable to any three-blade rotor 
at the same values of V, P, and H*. 

Four-blade rotor oscillatory aeroelastic derivatives, in stationary axes, 
contain harmonics of the product of number of blades and rotation rate b 0 

or 4p. Higher harmonics 8P, 12P are negligible and are not considered 

in this report. The theoretical variations of the 4P sine and cosine com- 
ponents of thrust coefficient divided by solidity, Cr p/a 5 derivatives with 
flap frequency P and advance ratio ^ are not shown. The derivatives with 

respect to 9 and 9_ , however, are much smaller and center about the 4P 
1c Is 

resonance of blade flapping instead of the 3P as was the case with the three- 
blade rotor. No thrust oscillation experimental data is available. 

It should be remembered that sines and cosines are relative to f = 0 
with blade number one pointed aft. 

The theoretical variations with P and |i of the 3P advancing rotating axes 
components of the 4P stationary axes hub moment aeroelastic derivatives with 
respect to 9^ and 9^ are shown in vector form for comparison with the equiv- 
alent components of the three-blade rotor, in Figures 83 and 84. Again the 
components are somewhat smaller and centered about different blade resonances. 
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Figure 84 . Advancing 3P Hub Moment Vectors Due to Unit Longitudinal 
Cyclic Pitch, at t = 0, for 4-Blade Rotors with V =5.0 



For comparison with experimental data, the pitch and roll components of 
the hub moment cyclic pitch derivatives are replotted versus advance ratio 
U* at the specific values of P tested and are shown in Figures 85 through 88. 
All data is for a blade Lock number V = 5.0. In addition, the hub moment 
component derivatives with respect to collective pitch 0 q and rotor angle- 
of-attack a are also shown in a similar manner and are compared with experi- 
mental data in Figures 89 through 92. 

Agreement of theory and experiment is quite good where the derivatives 
are of large magnitude. For the 7* 5-foot four-blade rotor, blade twist and 
precone were zero. Therefore rotor hub force residual harmonic components 
were theoretically zero. 
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Figure 85 . Hub Moment Oscillatory Aeroelastic Derivative With Respect to Lateral Cyclic Pitch, 3P 
Pitch Component at $ = 0 for 4-Blad.e Rotors With V = S.O 
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Figure 86. Hub Moment Oscillatory Aeroelastic Derivative With Respect to Lateral Cyclic Pitch 
3P Roll Component at = 0 for U-Blade Rotors With V = 5.0 
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Figure 87 . Hub Moment Oscillatory Aeroelastic Derivative With Respect to Longitudinal Cyclic 
Pitch, 3P Pitch Component at »j/ = 0 for 4-Blade Rotors with V = 5.0 
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Figure 88. Hub Moment Oscillatory Aeroelastic Derivative With Respect to Longitudinal Cyclic 
Pitch, 3P Roll Component at = 0 for 4-Blade Rotors With V =5.0 
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Figure 89 . Hub Moment Oscillatory Aeroelastic Derivative With Respect to Collective Pitch, 3P 
Pitch Component at if = 0 for 4-Blade Rotors With y = 5.0 
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Figure 90 . Hub Moment Oscillatory Aeroelastic Derivative With Respect to Collective Pitch 
3P Roll Component at = 0 for 4-Blad.e Rotors With V =5*0 












231 


Figure 92. Hub Moment Oscillatory Aeroelastic Derivative With Respect to Rotor Angle-of-Attack, 
3P Roll Component at ♦ = 0 for 4-Blade Rotors With V = 5.0 



Fixed shaft, free swashplate . - Fixed-shaft, free-swashplate tests at 
high advance ratio were performed only on the 33-foot 3-blade rotor-gyroscope 
system. For this reason analyses are limited to this rotor configuration — 
although the results are couched in a general fashion. 

In discussing free-swashplate gyroscope-stabilized systems, it becomes 
apparent that the rotation rate and size, or diametral inertia, of the gyro- 
scope and the feathering inertia of the blades as well as the blade forward 
sweep angle become important parameters. That is, as well as advance ratio u, 
lock number Y and flap frequency ratio P, the major parameters needed to de- 
scribe the fixed swashplate, fixed shaft response of the rotor, it is neces- 
sary to consider further nondimens ional parameters that control the free- 
swashplate behavior. 

The most notable component of the motion of such systems, especially 
those employing high speed gyroscopes, is the precesslve mode of the rotor - 
gyroscope combination. (The nutating mode natural frequency can be as high 
as POP and is therefore not significantly excited. ) 

The precessive mode advances and in stationary axes is of low frequency, 
typically of the order of P-1* If the ratio of its frequency to the shaft 
rotation rate can be preserved, say in a scaled model test in which the ad- 
vance ratio, flap frequency ratio, and Lock number of the blades have also 
been preserved then the nondimensional feathering response, control mode 
stability, and rotor loads of the full-scale rotor should be faithfully 
reproduced. 

The precessive or control mode natural frequency in stationary axes, 
neglecting the effects of blade flapping, is given in Reference 15 as: 
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where 
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control mode undamped natural frequency rad/sec 
rotor rotation rate rad/sec 
gyroscope rotation rate rad/sec 

2 

gyroscope diametral moment of inertia slugs ft 
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I + A L = blade feathering moment of inertia slugs ft 
o b 

mechanical advantage 
number of blades 

aerodynamic swashplate moment per unit swashplate tilt 
ft-lb/rad 

the ratio of control mode natural frequency to rotor rota- 
it is only necessary to preserve the ratio, 
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An approximate expression for f is given in Reference 14 for a three- 
blade rotor. It may be factored by ^ to make it approximately represent 
rotors with any numbers of blades. 
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Since ju is separately preserved, no account of it is required in the 
control mode preserving parameter. In addition, if the slope of the blade 
section lift curve is expressed as M a M and replaces the 2tt in the above ex- 
pression, then the part of "f" dealing with \i 
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may be dropped from the expression leaving the parameter which preserves the 
control mode frequency 
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f 
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Feathering Lock number 


It is also necessary to establish nondimens ional ratios for blade 
feathering friction and swashplate stationary axis damping C . These are 
as follows : 


and 
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C Ab 4 
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C q Ab n d 4 

b -- — per R 


where C = Swashplate feathering friction in rotating coordinates 
ft -lb/rad/ sec 

Cg = Swashplate damping in stationary axes, ft-lb/rad/sec 


In summary, the nondimensional motions, stability and response of any 
free -swashplate, fixed-shaft rotor are defined if the following nondimen- 
sional parameters describing the rotor-gyroscope system are known: 


Advance ratio 

Blade flap frequency ratio 
Blade Lock number 
Blade feathering Lock number 
Feathering friction coefficient 
Swashplate damping coefficient 
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P 

Y 
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It should also be noted that, although not necessary for stability, 
shaft-transmitted forces and higher harmonic blade motions and loads also 
depend on the number of blades b. 

The solution of the equations of motions of the free -swashp late, fixed- 
shaft system for motions and shaft transmitted forces follows: 

For the free-swashplate solution: 

■0» = e< = = 0 

<£ = $ = $ =0 

2 = Z = Z =0 

The swashplate tilt angles 0 and 4> are degrees of freedom, in contrast 
to the fixed-swashplate solution of the previous section. 

The two solutions differ mainly in the definition of the upper and lower 
portions of the vector. Since the swashplate tilt angles are now degrees 
of freedom, they are included in the {p^| vector: 
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Note that 
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The | vector of rotor loads will be moved to the left-hand side of 

equation (l6) with no change in sign, thereby becoming reaction moments. The 
vector of swashplate moments will be retained on the right-hand side, 
in which case they remain applied moments. Equation (l6) becomes: 
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(17) 


The preceding equation may be separated into two uncoupled parts. The 
lower portion becomes: 



(18) 


from which the displacement vectors may be determined. The known displace- 
ment vectors may then be used to get the response loads in the upper part of 
equation (17), as follows: 


W • • [ [ ir] M 



(19) 


The solution of equations (l8) and (19) proceeds in like manner to the 
prior section. The differences are minor and are due only to the differing 
matrix sizes and the somewhat different forcing functions of equations (6) 
and (18). 

Free-swashplate effects may be compared with those of the locked swash- 
plate only if the calculated cases are similar. It would not be reasonable, 
for example, to compare rotor shaft force derivatives, fixed and free, due 
to collective pitch or rotor angle -of -attack. They would be basically differ- 
ent due to the trimming effect of the free swashplate. 
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For a valid determination of the higher harmonic effects of the free 
swashplate on the shaft transmitted vibratory forces, the fixed-shaft free- 
swashplate system should be caused to change cyclic pitch, one component at a 
time, by judicial application of swashplate control moments. A comparison 
could then be made between the rates of change of shaft and swashplate 
forces with respect to unit changes in cyclic pitch for the fixed and free 
swashplates. Differences could then be logically ascribed to freeing the 
swashplate . 

The speculation that the swashplate motions, which should have 
occurred during the tests, were suppressed by high values of swashplate 
stationary axis damping was also investigated by calculating results with 

k = k = 0 as well as with the nominal values observed during bench 

, R. . s 
testing . 

Tip path motions and swashplate wobble: — The rates of change of non- 
dimensional tip path motions with respect to each of the cyclic pitch com- 
ponents were not measurably affected by freeing the swashplate, even in the 
absence of swashplate damping and feathering friction. In addition, the 
swashplate wobble was very small. 

This calculated result occurred at all conditions of rotor rpm and for- 
ward speed for which the system was stable. Inspections of the experimental 
data tended to confirm this theoretical finding. However, directly compar- 
able cases were difficult to find and the result was obscured by the free 
swashplate T s ability to neutralize random very low frequency perturbations in 
the flow field. 

The theory and experimental data indicated that mean and oscillatory 
aeroelastic force and moment derivatives for the stable rotor-gyroscope sys- 
tem were not significantly affected by freeing the swashplate. 

It should not be inferred from these results that no swashplate param- 
eter (feathering Lock number) and gyroscope inertia and rpm values exist 
which will reduce vibration. Preliminary studies by Dr. G. J. Sissingh 
(Reference 16) suggest that careful tuning of rotor -gyroscope systems can 
result in reduced vibration. No attempt, however, is made in this study to 
verify his predictions. 
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Free shaft, free swashplate. — Solution of the complete rotor-gyroscope - 
airframe equations for the motions, shaft -transmitted loads, and blade section 
forces is the objective of the free -shaft, free-swashplate calculations. 

These calculations are somewhat simpler than the previous analyses since 
none of the degrees of freedom are suppressed. 

The solution produces loads on free-flying aircraft in steady flight. 

The body degrees of freedom of the system are not restricted. 

This solution is the simplest of the three types of application. Again, 
it differs from its predecessors mainly in the definition of the j p} vector. 
In this instance, | , which is a vector of the locked-out degrees of 

freedom, does not exist. Thus: 



The solution will simplify to the form of equation (4) without the top 
row and first column, which is: 
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For the present applications, the vehicle will be controlled by moment 
inputs to the swashplate. Consequently there will be no loads applied to the 
hub such that the loads vector becomes: 
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0 




0 

0 

M e 

0 

0 

0 


Equation (19) is solved for the displacement vectors in like manner to 
equations (5) and (l8) of prior sections. 

No aeroelastic response derivatives may be determined for the free- 
flight cases. 
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SOLUTION OF IN -PLANE MOTIONS EQUATIONS 


Stability of In-Plane Motions 

Before undertaking the calculation of the steady oscillatory response of 
the in-plane blade -hub -body-spring system to its external forcing functions, 
it is advisable to examine its stability and modal behavior in the unforced 
state. This is done by arbitrarily letting the right-hand side (RHS) of the 
equation be zero. The left-hand set (LHS) or homogeneous equations may then 
be solved for its eigenvalues, or frequency and damping, and eigenvectors, or 
mode shapes. 

The damping characteristic of each mode of the system must be negative 
for stability and the system must be stable for realistic calculations of the 
oscillatory response to forcing functions. The natural or unforced frequen- 
cies of the modes vary with rotor rpm and since the modes are stable but 
lightly damped they experience resonance with rotor forces when their frequen- 
cies coincide with the ftb or 3 per revolution harmonics. The rpm's at which 
the in-plane modes intersect the 3P line are therefore important since they 
may lead to severe rotor loads. 

The following sections present the values of the rotor in-plane param- 
eters employed in the 33-foot 3-blade rotor analyses. Then the types of 
modes such a system possesses and how its frequencies vary with rpm is 
discussed. 


In-plane parameters for the 33-foot rotor. — 

M, 


Blade mass 

Blade c.g. radial position 
Position of blade pivot 


c.g. 


about 

shaft 

Ib shaft 

about 

its pivot 

^ b pivot 

nonrotating 

U) . 
IP. 


2.52 slugs 

9.00 ft 
2.50 ft 
253 Slugs ft 2 
156 slugs ft 2 
53»? rad/sec 
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Blade section drag coefficient 
Rotating hub mass 
Rotating hub moment of inertia 
Structural damping fraction 
Body effective lateral mass 
Body effective longitudinal mass 
Body lateral spring 
Body longitudinal spring 


Q o 

“r 

I R 

V 

M 

y 

M 

x 

K 

y 

K 

x 


= 0.011 

= 17.8 slugs 

= 6.72 slugs ft 

- 0 

= 66 slugs 

- 110 slugs 

- 68?, 000 Ib/ft 

- 235,000 lb/ft 


? 


In-plane modes and frequencies. — The number of modes equals the number 
of degrees of freedom since the body is attached to the wind tunnel by springs 
and is therefore not free. 


The modes of the system at zero rpm consist of one rigid body rotation of 
zero frequency, one collective in-plane mode of high frequency, in which the 
blades and hub rotate in opposite directions and four blade-body lateral and 
longitudinal modes. 

The four blade -body modes are of greatest interest and at zero rpm are 
real, or all elements in each mode keep the same relationship to one another 
throughout the period of oscillation. The elements change amplitude and sign 
but keep the same distribution. The four modes are sketched in Figure 93. 

The modes resemble two of those measured; however, the frequencies cal- 
culated are somewhat higher. 

Modes of the rotating rotor, in stationary axes, contain real and im- 
aginary parts. That is, the modes have different shapes at different times 
during the period of oscillation. This variation can best be seen by exam- 
ining the modal vectors in the imaginary plane. The actual shape at any 
instant of time is given by the projection of the vector components in the 
real axis. A construction of the real shape at regularly spaced intervals of 
time provides a "moving picture" of the mode throughout its period. 

The four rotor modes are shown in Figure 9h at rotor rpm of ? 50- Two 
modes are practically uncoupled body modes and two are the rotor advancing 
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f = 8.07 cps 


f = 11.16 cps 


BODY MODES 



y 



f = 4.55 cps 


f = 13.^5 cps 


ROTOR MODES 


Figure 94. In-plane Modes at ?50 RPM for the 33-Foot 3-Blade 
Rotor. Shaft Motion Exaggerated. 



and regressing modes. The fact that modes are advancing or regressing can be 
determined from the vectors in the imaginary (i-r) plane by noting that they 
keep constant magnitude while rotating counterclockwise or clockwise through 
the period of oscillation. The projection in the real axis of the advancing 
mode, for example, shows to be positive at the initial time so that the 

mode essentially appears as follows: 



y 


x 


One quarter of the period later ^ is positive and 
shape appears as follows : 


0 and the mode 



It is obviou... that the blades are precessing in the same direction as the 
rotor rotation or are advancing. 

The variation of the unforced rotor frequencies with rpm is shown in Fig- 
ure 95. It should be noted that modes Intersect the 3P line at rpm = 128, 

160, 220, and 275* Only the regressive mode intersection of 3P did not 
excite measurable oscillations in the rotor. Large chordwise oscillations of 
the blades were noted at the other intersecting values of rpm. 
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ROTOR ADVANCING MODE 



Figure 95* Unforced Body and Rotor Vibration Mode Frequency Variation 
with RPM for the 33-Foot 3-Blade Rotor 



BLADE LOADS 


The determination of blade section lift and flapping bending moment 
azimuthal distributions from calculated rotor external forces and response 
motions is discussed in this section. Theoretical flap bending distributions 
are shown compared to experiment at two radial stations on the 33-foot, 

3 -blade rotor. 

For the fixed-shaft 33-foot rotor both experimental data and theoretical 
results are analyzed so as to provide the rates of change of flap bending 
distribution with respect to cyclic pitch components and the residual bending 
distributions due to rotor fixed geometric parameters when cyclic pitch is 
zero . 

For the XH-51A rotor, experimental azimuthal distributions are avail- 
able from existing flight test conditions. These provided net or combined 
effects of cyclic pitch and rotor fixed geometry for comparison with theory. 
Correlation with the present theory was very poor and is not shown. 

Principles Involved in Blade Loads Determination 

Blade section lift and flap bending are found from the displacements, 
velocities, and accelerations of the degrees of freedom produced by the 
vertical motions differential equations when solved for steady oscillatory 
response to external forcing functions. 

Bending moments at blade stations are not found from the curvatures of 
the blade. Far too few modes are employed to make this a reasonable pro- 
cedure. Bending moments are found by integrating the aerodynamic and inertia 
section forces factored by the moment arms to the station from the station to 
the blade tip. The blade mode shape primarily affects inertia acceleration 
and centrifugal forces. The aerodynamic forces are largely independent 
of mode shape. For these reasons fairly accurate flapping bending moments 
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were predicted at blade stations near the root at high rpm and at all stations 
at low rpm despite the fact that only one flap bending mode was employed. 


The actual procedure employed to calculate the blade loads azimuthal 
distributions is indicated below for the three-blade rotor. A similar pro- 
cedure was used for the four-blade case. The blade-f orce-at- station-r 
column matrix is defined as follows: 


BM ' 
V 



FM 
k SL 


bending moment at station r 
shear at station r 
pitching moment 
feathering moment 
section lift at station r 


| F b} r = + [ G V r l] (n b | ♦ [V r - e) ] j\) 

+ [ i b <r) ](\}" [ < y r i* ) ] {*f} 

Blade forces are calculated from blade deflections, velocities, and 
accelerations and external forces. The matrices in the blade force equation 
contain the following terms 


Q T]( r 5 ’H 1 ’ ) 
QpOO) 


aerodynamic, displacement contribution to blade 
section forces 

aerodynamic velocity contribution to blade section 
forces 

aerodynamic and inertia external forces 
centrifugal contribution to section forces 
acceleration contribution to section forces 


Elements in the blade forces matrices are calculated in a way parallel 
to that employed in calculating coefficients in the differential equations. 

The vectors of blade b deflections ?! and fixed geometric shapes rj 

D 'F 

are defined as follows: 
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The blade "b" vector and its derivatives are determined from the overall 
rotor displacements, velocities, and accelerations relative to rotor coord- 
inates through the use of the restraining transformation matrices. 



and the motions in rotating coordinates depend on the motions in stationary 
axes through the time dependent transformation [ t] . 

b r = T e s 

e R - TB S ♦ TB S 

e R = T 0 S + 2T B s + T B s 

and 8 is the vector of degrees of freedom of the equation in stationary 

o 

axes; the output of the steady oscillatory calculations. 



and is itself a function of azimuth. 

B = 8 + cos 3^ + 8_ sin 3^ + B/- cos 6+ + sin 6^ 

o 3c 3s be os 
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Blade Flap Bending Moments, 33 -Foot Rotor 

On the 33 -foot 3 -blade rotor, flap bending moments were measured at 

P = .217 and .596 (stations m 3 and ll8 inches respectively) at 80 knots 

forward speed over a range of rpm. At ^ = . 217 , moments on the three blades 

were analyzed and presented in nondimensional form. The measurements showed 

a common basic behavior and provided a check on each other despite fairly 

r 

large differences between the blades. At — = .596, bending moment measure- 
ments were available only on blade number 2. 

The test cases analyzed are shown in the following table. 


Forward 

Speed 
( Knots '■ 

RPM 

M= 

P 

v 

80.59 

197.2 

.399 

1.35 

4.57 

82.76 

153.2 

.528 

1.53 

4.57 

82.78 

100.2 

co 

O 

00 

2.03 

4.57 

82.68 

72.1 

1.121 

2.64 

4.57 

82.88 

38.0 

2.132 

4.67 

4.57 


At each radial station for each test condition, bending moment distri- 
butions caused by approximately a dozen combinations of cyclic pitch were 
analyzed. Best fit planes (rms fit) of each harmonic component versus the 
two cyclic pitch angles supplied rates of change of bending moment harmonic 
component with respect to each cyclic pitch component and a residual value 


From these, bending moment azimuthal distributions per unit value of 

each cyclic pitch component and the residual were prepared. These were then 

compared with theoretical values. All bending moments were presented in the 

form - Ch w where C. is defined: 
a D.m. b.m. 

c __ bending moment 

b ' m * p(f)R ) 2 hr 2 r 
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For puiposes of comparison the azimuthal distributions of bending moment at 
the two radial stations due to longitudinal cyclic pitch G^ s are shown in Fig- 
ures 96 through 105 for the five tested cases in the order given in the table. 

It is interesting to note, at high rpm, in the first two cases, that the 
bending moment at ^ = -596 is not correctly predicted at all, where, as the rpm 
reduces over the last three cases, the agreement between theory and experiment 
improves. This is thought to be due to the single parabolic mode shape employed 
to represent blade flapping deflection. The tip region, at high rpm, would the- 
oretically be deflected upward too far and would induce a negative bending moment 
at a point in the azimuth where a positive bending moment occurred physically. 

At low rpm, on the other hand, tip deflections and centrifugal forces are 
very small and most of the blade bending moment is produced by aerodynamics. 

Figures 96 through 105 also show the transition of the distribution from an 
almost sinusoidal form at P = 1.35 and 1.53 to heavy two-per-revolution oscilla- 
tions at P = 2.03. At the latter condition the blade is in damped resonance with 
its 2P exciting forces. At the P = 2.64 condition light three per revolution 
oscillations become evident and at P = 4.67 four per revolution become evident. 

Figures 106 through 115 show the variation of the azimuthal distribution pro- 
bending moment at — = .217 and .596 due to unit lateral cyclic pitch 0 over the 
n -LC 

range of tested cases. Remarks regarding the variations with longitudinal cyclic 
pitch apply in general. However, the change in phase due to the change in cyclic 
pitch component may be noted. 

Figures 11 6 through 125 show the bending moment azimuthal distribution pro- 
duced by precone (3 =2.25 degrees, blade twist 0,R = -9.^-3 degrees, and collec- 

o t 

tive pitch 0 „ =1.5 degrees, with cyclic pitch and angle of attack zero = 

• f 5-R -tc 

0 = a = 0. These plots include the centrifugaUy induced bending moment due 

J. s 

to precone and collective blade flapping, , and at high rpm, the first 
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two cases, are suspected of causing the deviation of the mean bending moment 
from the experimental values. At low rpm the effect is not quite as 
pronounced . 

The residual flap bending moment distribution roughly resembles the 
negative of the bending moment due to 0 . Trimming bending moments to zero 
therefore could be approximately accomplished by an application of positive 



Reference 3? Figure 5 show + 6 to be the larger of the cyclic pitch 

components required to trim hub moment to zero, although fairly large values 

of - 0_ are also required. 

Is 

In summary, the theory predicted the essential features of the bendinr 
moment distributions at the inboard station at the five combinations of F 
and p tested at a Lock number V = 4.57. 
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Figure 96 . Blade Flap Bending Moment per Unit ©is, 33 -Foot Rotor 
H- = -399, P = 1-35, = 4.57. r/R = .217 
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Figure 97- Blade Flap Bending Moments per Unit 0p s , 33-Foot Rotor, 
H- = -399, P = 1-35, = 4.57. r/R = .596 
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Figure 100, Blade Flap Bending Moment per Unit Bis? 33-Foot Rotor 
Hi = .808, P = 2.03, V = 4.57. r/R = .217 
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Figure 101. Blade Flap Bending Moment per Unit 0i S3 33-Foot Rotor 
P - .808, P = 2.03, Y = 4.57. r/R = .596 
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Figure 102. Blade Flap Bending Moment per Unit 0p s , 33-Foot Rotor 
M- = 1.121, P = 2.64, Y = 4.57. r/R = .217 
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Figure 103. Blade Flap Bending Moment per Unit 0i s , 33 -Foot Rotor 
4 = 1.121, P = 2.64, Y = 4.57. r/R = .596 
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Figure 1CU. Blade Flap Bending Moment per Unit G]_ S5 33-Foot Rotor 
ut = 2.132, P = 4.67, Y = 4.57 r/R = .217 
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Figure 105. Blade Flap Bending Moment per Unit Gp s , 33-Foot Rotor 
n = 2.132, P = 4.67, Y = 4.57 r/R = .596 
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Figure 108. Blade Flap Bending Moments Per Unit 33-Foot Rotor 

M. = .528, P = 1.53, Y = 4.57. r/R = .217 
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Figure 109. Blade Flap Bending Moments Per Unit Sic, 33-Foot Rotor 
H = .528, P = 1.53, Y = 4.57- r/R = .596 
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Figure 110. Blade Flap Bending Moments Per Unit 33-Foot Rotor 

M- = .808, P = 2.03, Y = 4.57. r/R = . 211 
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Figure Him Blade Flap Bending Moments Per Unit ®ic, 33-Foot Rotor 
M> = .808, P = 2.03, Y = 4.57. r/R = .596 
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Figure 120. Blade Residual Flap Bending Moments, 33-Foot Rotor, 
H. = .808, P = 2.03, Y = 4.57- r/R = .217 
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Figure 121. Blade Residual Flap Bending Moments, 33-Foot Rotor 
M> = .808, P = 2.03, Y = 4.57. r/R - .596 













CONCLUDING REMARKS 


1. A system of linear ordinary differential equations with periodi- 
cally varying coefficients, forced by functions containing oscillation of 
the same period, has been derived to represent the flapping and associated 
motions of a hingeless rot -r, free-swashplate, f ree-a i rframe configuration 
operating at advance ratio greater than .3 (u > ♦3) and with its ratio of 
blade flap frequency to rotor rotation rate greater than 1.3 (P > 1.3) • 

2. The system of equations has been solved for its steady oscillatory 
response to the periodic forcing produced by cyclic pitch, collective pitch 
and angle-of-attack, with some of its degrees of freedom suppressed. The 
shaft-transmitted moments obtained from the response were resolved into 
mean and oscillatory aeroelastic derivatives and approximately agreed with 
experimental data reduced and analyzed so as to provide equivalent 
derivatives . 

3- Shaft moment oscillatory aeroelastic derivatives agreed best with 
experimental data when the blade natural flap frequency was in the vicinity 
of the two-per-revolution excitation, in rotating axes, for the three-blade 
rotor and in the vicinity of the three-per-revolution excitation for the 
four-blade rotor. 

4. Omission of the harmonic components of the coefficients in the 
differential equation of vertical motion yielded approximately correct 
mean aeroelastic derivatives, at all values of flap frequency ratio pro- 
vided the advance ratio was less than about .8 (p < .8). 

At advance ratio greater than .8 omitting the periodic components of 
the coefficients caused the mean aeroelastic derivatives to be greatly in 
error. 

5. In addition to shaft forces, the equations of motion have been 
solved for the azimuthal distribution of blade flapping bending moment. 

The bending moment distributives were computed for unit values of cyclic 


268 



pitch and compared with experimental distribution deduced from measured 
data in an equivalent form. 

Despite the fact that blade flap bending moments are sensitive to the 
shape of deflection modes employed in the theory, and the theory of this 
report used only a single parabolic deflection degree of freedom, bending 
moment distributions calculated were quite similar to those deduced from 
the experimental data. 

Blade root region flap bending was approximated correctly at all tested 
values of blade flap frequency ratio (1.3 < P < 5) at advance ratio greater 
than .5 (u > .5) • Blade outer region bending moments, however, became 
realistic only at large values of flap frequency ratio (P >3)* 

6. A set of linear ordinary differential equations describing the 
in-plane motion of hingeless rotor blade-hub-body-spring system were de- 
rived and applied to the 33“foot 3-blade rotor. Theoretical modal char- 
acteristics were compared with experimental. 

Means of determining the in-plane forcing functions based on the 
motions and forces of the vertical equations were discussed. 
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